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I. INTRODUCTION: HOW TO 
CHARACTERIZE ENTANGLEMENT? 

Inseparability of quantum states of composite systems, 
discovered in the early days of quantum mechanics by A. 
Einstein, B. Podolsky, and N. Rosen [lj and named "en- 
tanglement" ( "Verschrankung" ) by E. Schrodinger 0, 
became one of the central concepts of contemporary 
physics during the last decade. Entanglement plays now 
a vital role within quantum information science [3j , repre- 
senting both the defining resource for quantum commu- 
nication - where it enables, in particular, non-classical 
protocols such as quantum teleportation 4] and it leads 
to enhanced security in cryptographic tasks Q - and a 
key ingredient for determining the efficiency of quantum 
algorithms and quantum computation schemes [g, 0] • ^ n 
addition, studies of entanglement have also proved be 
relevant to fields as different as atomic physics Q , quan- 
tum chaos |9| . fllL quantum phase transitions 
Q IS H ll6lTl7l TTa TTa l20| . and quantum networks 

According to the original definition, the description of 
entanglement relies on a specific partition of the compos- 
ite physical system under consideration. However, such 
a system can often be decomposed into a number of sub- 
systems in many different ways, each of the subsystems 
possibly being a composite system by itself. In order to 
avoid ambiguity, given a partition of the composite sys- 
tem into n subsystems, we call each of them an "element" 
and characterize it by a single, possibly collective, quan- 
tum number. Thus, the composite system is a collection 
of the elements. We call this collection an "assembly" 
in order to avoid confusion with an "ensemble" , which 
is usually understood as a set of all possible realizations 
of a many-body system with an associated probability 
distribution over these realizations. The i-th element is 
assumed to have Hilbert space of dimension di. Qubit, 
qutrit, and qudit are widely used names for two-, three-, 
and d-level elements (with di = 2, di = 3, and di = d, 
respectively). 



The term "entanglement" has a transparent qualita- 
tive meaning: A pure state of an assembly is entangled 
with respect to a chosen partition when its state vector 
cannot be represented as a direct product of state vec- 
tors of the elements. This notion can also be extended to 
generic mixed quantum states, whereby entanglement is 
defined by the inability to express the assembly density 
operator as a probabilistic combination of direct prod- 
ucts of the density operators of the elements. Intuitively, 
one expects that in the presence of "maximum" entan- 
glement the states of all subsystems comprising the 
overall system are completely correlated in such a way 
that a measurement performed on one part determines 
the states of all other parts. 

The question is: how to quantitatively characterize 
entanglement for an assembly of many elements? One 
would like to have a measure ranging from zero for 
the product state to the maximum value for the maxi- 
mally entangled state. This can be easily accomplished 
for the bipartite setting n — 2 that is, for an assem- 
bly consisting just of two distinguishable elements A 
and B, each of arbitrary dimension. In this case, the 
von Neumann Tr[p^ log pa] or linear Tr[p^] entropies 
based on the reduced density operator of either element, 
e.g., pa = Ttb[\^) (^l], may be chosen as entanglement 
measures for a pure state \^} of the assembly. How- 
ever, already for the tripartite case n = 3, characteriz- 
ing and quantifying entanglement becomes much harder. 
In fact, there are not one, but many different charac- 
teristics of entanglement and, apart from the question 
"How much?" , one has also to answer the question "In 
which way" |23j, |24j are different elements entangled? In 
other words, apart from separability criteria |25|. one 
also has to introduce inequiyalent entanglement mea- 
suresJMJplJp, |29l |30l l3ll lp| and entanglement classes 

[iiiifiMifr 

The identification of appropriate measures is not 
unique, and is mostly dictated by convenience. The 
choice of classes may have a more solid ground based on 
group theory, since their definition is related to groups of 
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local operations. These are operations applied individu- 
ally to each of the elements and forming a subgroup of 
all possible transformations of the assembly state vector. 
One may resort to group-theoretical methods, which al- 
low one to construct a set of invariants |33l 13(1 1371 l3q 
under such local operations. Under the action of local 
transformations, which can be either unitary or, in the 
most general case, simply invertible, the state vector of 
the system undergoes changes, still remaining within a 
subset O (an orbit) of the overall Hubert space Ti. The 
dimension of the coset H./O, that is the number of in- 
dependent invariants identifying the orbit, can be easily 
found for a generic quantum state. Still, there are sin- 
gular classes of orbits that require special consideration. 
Their description depends critically on the number of el- 
ements and on the detailed structure of the assembly. 

The numerical values of the complete set of invariants 
may be chosen as the "orbit markers" that provide one 
with an entanglement classification, although the choice 
of this set is not, in general, unique. It turns out that gen- 
erally accepted measures of entanglement, such as con- 
currence [23 for two qubits (n = 2, di = 2), and 3-tangle 
[23| for three qubits (n = 3, di = 2), are such invari- 
ants [3^, • Generalization of these measures |3^| to 
n = 4, di = 2 provides a connection between measures 
and invariants characterizing different classes. Their con- 
struction is not an easy task, and one needs to identify 
the invariants that are able to distinguish inequivalent 
types of entanglement |4*flj . Moreover, the invariants are 
usually high-order polynomials of the state amplitudes, 
with the maximum power growing linearly with the num- 
ber of elements in the assembly. Therefore they rapidly 
become rather awkward |38l l4f| . A partial albeit not un- 
ambiguous connections between the measures and classes 
are established by the requirement jl^ f° r a measure to 
behave as a so-called monotone. This means it should 
be non-increasing on average under the action of non- 
unitary invertible transformations of the elements, also 
known as the family of Local Operations and Classical 
Communication, LOCC |43|. The complete classification 
problem remains unsolved even for relatively small as- 
semblies (see e.g. 01 f° r recent results on the the five- 
qubit system). 

Apart from polynomial-invariant constructions, other 
schemes have been proposed to describe multipartite en- 
tanglement, including those based on generalization of 
Schmidt decomposition [34I EBl l46l Efl. on invariant 
eigenvalues |33|, on hyperdeterminants l35[ . and on ex- 
pectation values of anti- linear operators |4Cj |. However, 
none of these suggestion has been fully tested for n > 4 
qubits or more than three qutrits (n > 3, d = 3) |48[ . 
Moreover, for the orbits of general invertible local trans- 
formations, the complete sets of invariants are unknown 
|49j for assemblies consisting of n qudits (di = d) if n > 3 
and/or d > 4. Still, a number of physically reasonable 
suggestions |28|, H2, Om |5j| for entanglement characteri- 
zation have been attempted. 

In this paper, we propose a different approach to en- 



tanglement characterization. We focus exclusively on the 
case where an assembly in a pure quantum state con- 
sists of distinguishable elements, leaving generalizations 
to mixed states and to indistinguishable elements for fu- 
ture studies. Our main aim is to construct extensive 
characteristics of entanglement. Thermodynamic poten- 
tials linearly scaling with the number of particles in the 
system offer examples of extensive characteristics widely 
employed in statistical physics. The free energy given by 
the logarithm of the partition function is a specific impor- 
tant example. We will introduce similar characteristics 
for entangled states in such a way that their values for a 
product state coincide with the sum of the corresponding 
values for unentangled groups of elements. This tech- 
nique is based on the notion of nilpotent variables and 
functions of these variables. An algebraic variable x is 
called nilpotent if an integer n exists such that x n = 0. 
In our case, these variables are naturally provided by 
creation operators, whence the logarithm function trans- 
forming products into sums plays the central role in the 
construction. 

Our approach is based on three main ideas: (i) We 
express the state vector of the assembly in terms of a 
polynomial of creation operators for elements applied to 
a reference product state, (ii) Rather than working with 
the polynomial of nilpotent variables describing the state, 
we consider its logarithm, which is also a nilpotent poly- 
nomial. Due to the important role that this quantity 
will play throughout the development, we call this quan- 
tity the nilpotential henceforth, by analogy to thermody- 
namic potentials. (Hi) The nilpotential is not invariant 
under local transformations, being different in general for 
different states in the same orbit. We therefore specify 
a canonic form of the nilpotential to which it can be 
reduced by means of local transformations. The nilpo- 
tential in canonic form is uniquely defined and contains 
complete information about the entanglement in the as- 
sembly We therefore call this quantity the tanglemeter. 
The latter is, by construction, extremely convenient as 
an extensive orbit marker: the tanglemeter for a system 
consists of several not interacting, unentangled groups of 
elements equals the sum of tanglemeters of these groups. 

Let us briefly explain these ideas, in the simplest ex- 
ample of n qubits, which will be discussed in detail in 
Sect. II. An assembly of n qubits is subject to the 
sn(2)i®. . .ffisu(2)i . . .©su(2)„ Lie algebra of local trans- 
formations [52^. As a reference state, we choose the Fock 
vacuum that is, the state |0) = |0, 0, . . . , 0) with all the 
qubits being in the ground state. An arbitrary state of 
the assembly may be generated via the action of a poly- 
nomial F(a~l) in the nilpotent operators af on the Fock 
vacuum. Here, the subscript i enumerates the qubits, 
and the operator creates the state |1) out of the state 
|0). Evidently, (cr^) 2 = 0, since the same quantum state 
cannot be created twice. The family of all polynomials 
F(o~^~) forms a ring. We note that anticommuting nilpo- 
tent (Grassmann) variables are widely employed in quan- 
tum field theory [5j| and in condensed matter physics 
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|54j . However, the nilpotent variables introduced here 
commute with one another. 

In order to uniquely characterize entanglement, one 
must first select a convenient orbit marker. Following the 
idea of Ref. ^5|, we take as such a state |^ c )c) ly m S m 
the orbit O, which is the "closest" to the reference state 
|0) in the inner product sense, that is |(0 |\I/ C ) Q | = max. 
Once the state |^ / c )o i s found, it is convenient to impose 
a non-standard normalization condition |(0 |vt , c )ol = 1- 
We call the resulting state I^^q canonic. The canonic 
state is associated with the canonic form of the polyno- 
mial F c (af), which begins with a constant term equal to 
1. 

We mainly work not with F c by itself, but with the 
tanglemctcr a nilpotent polynomial f c = In F c that can 
be explicitly evaluated by casting the logarithm function 
in a Taylor series of the nilpotent combination F c — 1. 
Since (c^~) 2 — 0, this series is a polynomial containing 
at most 2™ terms. Both the tanglemeter and nilpotential 
(/ = In F) resemble to the eikonal, which is the logarithm 
of the regular semi-classical wave function in the position 
representation, multiplied by —i. The difference is that 
in our case no approximation is made: / represents the 
logarithm of the exact state vector. 

The nilpotential / and the tanglemeter f c have sev- 
eral remarkable properties: (i) the tanglemeter pro- 
vides a unique and extensive characterization of entangle- 
ment; (ii) a straightforward entanglement criterion can 
be stated in terms of the cross derivatives d 2 f / 'daf da^ '; 
(Hi) the dynamic equation of motion for / can be writ- 
ten explicitly and, suggestively, in the rather general case 
has the same form as the well-known classical Hamilton- 
Jacobi equation for the eikonal. 

The paper is organized as follows. In Sect. II, we an- 
alyze in detail an assembly of n qubits in terms of the 
nilpotent polynomials F and /. We extend the notion 
of canonic forms to the group of reversible local transfor- 
mations SL(2, C) and introduce the idea of entanglement 
classes. We conclude the section by presenting expres- 
sions relating the coefficients of Fc and fc with known 
measures of entanglement. To avoid confusion, we note 
that the subscript c corresponds to su-canonic forms in 
contrast to C, which corresponds to sZ-canonic forms. 
Details of the calculations and some proofs are given in 
Appendices A and B along with graphic representations 
of the entanglement topology. 

In Sect. Ill, we consider the evolution of the nilpotent 
polynomials under the action of single-qubit and two- 
qubit Hamiltonians, and derive an equation of motion for 
the nilpotential, which is distinct from the Schrodinger 



equation. For one important particular case able to sup- 
port universal quantum computation |55j . we show that 
this equation has a form of the classical n-dimensional 
Hamilton- Jacobi equation. Describing quantum dynam- 
ics in terms of nilpotentials suggests a computational al- 
gorithm for evaluating tanglemeters, which can be per- 
formed by dynamically reducing the polynomials to the 
forms canonic under either SU (2) or SL(2, C) local trans- 
formations. In the example of a four-qubit assembly, we 
explicitly illustrate how to identify the resulting entangle- 
ment classes. This technique yields entanglement classes 
consistent with the results of Ref. [24[ ■ The explicit anal- 
ysis of these classes as well as details of derivation of the 
equation of motion are given in Appendices B and C. 

In Sect. IV, we extend our technique to assemblies 
of ci-level elements - starting from qutrits [56L |5?J ■ The 
Cartan-Weyl decomposition of the su(d) algebras sug- 
gests a natural choice of nilpotent variables for qudits. 
For each element, we have d — 1 = r variables represent- 
ing commuting root vectors from the corresponding Lie 
algebra, which has rank r. For the illustrative case of 
two and three qutrits, we discuss possible choices of the 
canonic forms of the nilpotent polynomials. We further 
extend the approach to the case where the assembly par- 
tition may change as a result of the merging of elements, 
such that the new assembly consists of fewer number of 
elements with di ^ dj, and consider transformations of 
nilpotent polynomials associated with such a change. Fi- 
nally, we address a situation encountered in the frame- 
work of generalized entanglement |58l 59] , where the rank 
r of the algebras of allowed local transformations is less 
than d — 1. In other words, while the assembly is still 
assumed to be composed of a number of distinct ele- 
ments, the group of local operations need not involve 
all possible local transformations. In such a situation, 
the proper nilpotent variables are more complicated than 
a + . In particular, they may have non- vanishing squares 
etc, with only di powers vanishing. In addition, unlike 
in the conventional setting, entanglement relative to the 
physical observables may exist not only among different 
subsystems, but also within a single element. 

We conclude by summarizing our results and discussing 
possible developments and future applications of nilpo- 
tent polynomials and the tanglemeter. 

II. ENTANGLEMENT CHARACTERIZATION 
VIA NILPOTENT POLYNOMIALS 

Consider n qubits in a generic pure state |^), 



l*> = 2^ '0fe„fc„_ 1 ...fci \K, K-x, ■ ■ ■ i h) = ipoa...o |o,o, . . .0) 

{i,}=0,l 

+ V10...0 |1, 0, ... 0) + Vol. .o |0, 1, . . . 0) + . . . + 11, 1, • • • 1) 



(1) 
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specified by 2 n complex amplitudes tpk„k n -i-.kii i- e - by 
2«+i rea j numbers. The index fcj = 0, 1 corresponds to 
the ground or excited state of the i-th qubit, respectively. 
When we take normalization into account and disregard 
the global phase, there are 2" +1 — 2 real parameters char- 
acterizing the assembly state. 

It is natural to expect that any measure characteriz- 
ing the intrinsic entanglement in the assembly state re- 
mains invariant under unitary transformations changing 
the state of each qubit. A generic SU(2) transformation 
is the exponential of an element of the su(2) algebra, 

U = exp[i(a x P x + a y P y + a z P z )\ , (2) 

where a x , a y , and a z are Pauli matrices. It depends on 
the three real parameters P x , P v , and P z . Such a trans- 
formation changes the amplitudes ^fc n fc n _i...fci i n Eq. QJ, 
but preserves some combinations of these amplitudes - 
the invariants of local transformations. Thus, local trans- 
formations move the state along an orbit O, while the 
values of the invariants serve as markers of this orbit. 

The first relevant question is: What is the maximum 
number of real invariants required for the orbit identifica- 
tion, hence, for entanglement characterization? A generic 
SU(2) transformation represented by Eq. @ depends on 
three real parameters. Therefore, for n qubits the dimen- 
sion of the coset Tt/O, that is the number of different 
real parameters invariant under local unitary transfor- 
mations, reads HEE3 

D su = 2" +1 - in - 2 . (3) 

Mathematically, the counting in Eq. © corresponds to 
the number of the invariants of the group ®iSUi(2) ® 
U(l), where the factor £7(1) describes the multiplication 
of the wave function on the common phase factor and the 
usual normalization condition (vP|vf) = 1 for the wave 
function is imposed. It is more convenient for us to re- 
formulate the same problem as seeking for the invariants 
of ®iSUi(2)®C* , where C* is the group of multiplication 
by an arbitrary nonzero complex number and no normal- 
ization condition on the wave function is imposed. This 
will allow us to choose a representative on the orbit of O 
with nonstandard normalization choice (^c|0) = 1. 

To be precise, the counting Eq. is true for n > 2 
while the case n = 2 is special: in spite of the fact that 
2 3 — 3 • 2 — 2 = 0, there is a nontrivial invariant of local 
transformations for two qubits. It has the form 

I = TpOQlpll ~ IpOllpW- (4) 

For a three-qubit system, five independent local invari- 
ants exist |39j, namely three real numbers 

h = Tpijkr ijP ^mn P r mnk , (5) 

and the real and the imaginary part of a complex number, 
h + ih = ^fcV^VWpV-" 1 ^ (6) 



Here, ■0 yfe = e M e-" e kk ipi'j'k', with the summation over 
repeated indexes taking values and 1 implicit, ip* l: > k 
denotes the complex conjugate of ipijk, and e !! is the 
antisymmetric tensor of rank 2. The quantity 2|/4 + 
is also known by the name residual entanglement or 3- 
tangle r Refs. |2<lls3- 

Similar invariants can still be found for a four-qubit 
system. However, with increasing n, the explicit form of 
the invariants becomes less and less tractable and con- 
venient for practical use. Moreover, no explicit physi- 
cal meaning can be attributed to such invariants. We 
therefore suggest an alternative way to characterize en- 
tanglement, which is based on: (i) Specifying the canonic 
form of the state that unambiguously marks an orbit; (ii) 
Characterizing this state with the help of coefficients of 
a nilpotent polynomial; (Hi) Considering the logarithm 
of this polynomial, the tanglemeter. Thus, we construct 
extensive invariants of local transformations as the co- 
efficients of the tanglemeter that is, nilpotential of the 
canonic state. 

In this section, we proceed with illustrating the main 
technical advantages of our description within the qubit 
setting, deriving an entanglement criterion, and explain- 
ing how the invariants constructed by our method are 
related to existing entanglement measures. We also an- 
alyze a case important for certain applications involving 
indirect measurements, where it is natural to consider a 
broader class of local transformations constrained only 
by the requirement of unit determinant. Specifically, we 
focus on the set of stochastic local operations assisted by 
classical communication 43], which is widely employed 
in quantum communication studies and protocols. For 
qubits, such transformations are known as SLOCC maps 
|3lj . These operations do not necessarily preserve the 
normalization of state vectors. However, as suggested 
by Theorem 1 of Ref. [3j|, after a proper renormaliza- 
tion they are described by the complexification sl(2, C) of 
the su(2) algebra, such that the parameters (P x ,Pf , P z ) 
specifying the transformation of Eq. J2J) on each qubit are 
now complex numbers. The corresponding real positive 
invariants of local SX(2,C) transformations are mono- 
tones. 



A. Canonic form of entangled states 

In order to unambiguously attribute a marker to each 
orbit, we specify a canonic form of an entangled as- 
sembly state. To this end, we first identify a refer- 
ence state |0) as a direct product of certain one-qubit 
states. The latter can be chosen in an arbitrary way, 
but the choice |0), with the lowest energy level occu- 
pied, is the most convenient to our scope. Thus, the 
reference state reads |0) = |0, ...0). Drawing paral- 
lels with quantum field theories and spin systems, we 
will call |0) the "ground" or "vacuum" state. Then, 
following a suggestion of Ref. 45], by applying local 
unitary operations to a generic quantum state |\&), we 
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can bring it into the "canonic form" \$f c ) correspond- 
ing to the maximum possible population of the reference 
state |0). In other words, we apply a direct product 
Ui ® • • • <8> U n of transformations as in Eq. @ to the state 
vector |\ff), and choose real parameters (Pf , Pf , Pf) to 
maximize |(0| U\ ® . . . (g> U n |^)| 2 - The transformation 
t/i ® . . . ® f/„ satisfying this requirement can be seen to be 



unique up to phase factors multiplying the upper states 
of each qubit. Modulo this uncertainty, the canonic state 
\^c) = f/i®. • -®U n l^) can serve as a valid orbit marker. 

In fact, a generic unitary transformation of the ru- 
th qubit, chosen in the form exp[ia z <p m ] exp^tr^gm] 
exp[ier z ifim], 4>m, grmfm € K, equivalent to Eq. (J2J, re- 
sults in 



ip...k m =o... — > ^...fe m =0...e 10m IVm cosg m + i'0...fc„ l =i...e I¥>m I0m sin5 m , 

</>... fcm =i... -» ^...fc m =i...e i0m+iVm cos ffm + iV... fcm =o...e i0 "'- iv ™ sin ffm . (7) 

I 

Let us consider a generic infinitesimal local transforma- Eq. (J7J in series in g m up to the second order, one obtains 
tion of the state in the canonic form. By expanding 



O 



% L ) 1 /'o...o + iX!-9™ e2lV '"^o...fc m =i...o-X! 9mgie 21 



fe m = l...fc; = 1...0 



r 

(8) 



for the amplitude of the ground state. Since the parame- 
ters of the transformation are arbitrary, the condition of 
maximum ground-state population |(0|<I')| 2 implies that 
the linear term in Eq. © vanishes, 

ipo...k m =i...o = , Vm. (9) 

This gives n complex conditions and implicitly specifics 
2n out of in real parameters of the local transforma- 
tion that maps a generic state to the canonic form. The 
remaining n parameters may be identified in a generic 



case with the phase factors e -1 ^ r+ * 3r ), where one can set 
cj) r = without loss of generality. Two remarks are in 
order. 

(i) Special families of states of measure zero in the as- 
sembly Hilbert space may exist for which the system of 
equations @ is degenerate and specifies less than 2n pa- 
rameters. The simplest example for n — 2 is the Bell 
state, with ?/>oo = ipn — l/v2 and Voi = V'io = 0. The 
combination of two transformations of the form gives 
a state with the amplitudes 



#1 



V2 

e i(tfe-0i) 
e i(0i+0 2 ) 



V2 



e '(^i+v 2 ) singi cos (72 + e 1 ^ 1 ^ 2 -* cos sin (72 



(10) 



One can see that the conditions 7/^ = and ip[ — 0, y>i + (^2 = t/2 with arbitrary ^1,2)- 

are not independent: They both give 01 + 02 =0 and r™ u ., t ,,. . , , , , , 

_ •; 1 he orbit 01 this special state has 4 parameters. On the 

<Pi+<P2 = with arbitrary 6\ 2 (or, equivalently, q\ — 02 = ,1 v. j f • a j. 1 1 1 

r r ■ ' ^ ax other hand, for a generic canonic state with Vii/Voo — a i 
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| a | 7^ 1, the conditions ip' 01 = ip' 10 = imply g L = g 2 = 0, 
and only two phases (^i^ + Vi^ are arbitrary, whereas the 
transformed state being independent of the differences 
4>i,2 ~ 921,2 m this case. 

When n grows, the pattern of such special classes of 
states becomes more and more complicated. These fam- 
ilies resemble "catastrophe manifolds" where infinitesi- 
mal variation of the state amplitudes tp result in a finite 
change of the local transformations reducing the state to 
the canonic form. Here, we shall not discuss this further 
and restrict ourselves to the generic case. 
(ii) As noticed in Ref. 39], the conditions © are nec- 
essary but not sufficient in general for the state to have 
maximum ground-state population |(0|^)| . For exam- 
ple, an n = 2 state with |^n| > IV'ool and V'oi = i'w = 
does not have the maximum ground-state population, 
although it satisfies Eq. 10: when \ipn\ starts to ex- 
ceed IV'ooIj finite "spin-flip" operations must be applied 
to both qubits to reduce the state to the canonic form. 

For a generic n-qubit assembly state, the canonic form 
is unique up to n phase factors <p m + ip m , and the 
state may be characterized by 2" — n — 1 complex ra- 
tios a fc „fc„_ 1 ... fcl = ipk n k n - 1 ...ki/'>Po...o with J2 m k ™ > 
whereas the amplitude of the vacuum state ipo...o, after 
being factored out, specifies the global phase and the 
normalization. We disregard these factors and normalize 
the wave function such that the amplitude of the ref- 
erence state ipo.,.0 is set to unity. Then the parameters 
c*fc„fc„„i...fci correspond to the amplitudes of the assembly 
states where at least two elements are excited. The num- 
ber of real parameters characterizing the canonic form 
equals 2 n+1 — 2n — 2. It is worth mentioning that all 
|ctfc„fc„_i...fci | are invariant and, moreover, in the case n > 
3, the ratios ak n k n -i..M a l n l n -i..M/ a k! n k' n _ 1 ...k' x a l>J'^_ 1 ...l' x 
are invariant if for each m one of two conditions k' m = k m , 
l' m = l m or l' m = k m , k' m = l m are satisfied. By specifying 
n factors, we arrive at the bound Eq. J3J for the maxi- 
mum number of invariants characterizing entanglement. 

Indeed, by an appropriate choice of the phase factors 
in Eq. Q, one can make a set of n non-zero amplitudes 
Q!fc„fc n _i...fei real and positive. For example, for a generic 
orbit one can make n amplitudes otk n k n -i---k x correspond- 
ing to the next-to- highest excited states real and positive, 
with J2 m k m = n—1. In Fig.^i) we illustrate this for the 
simplest case n = 3 where the coefficients aon, «on, and 
aon are chosen to be real and positive. As mentioned, 
the case n = 2 is special, since the next-to- highest excited 
state amplitudes coincide with the first excited ones that 
vanish due to the requirement of Eq. @. Thus, a single 
parameter an characterizing entanglement can always 
be chosen real and positive, in accordance with Eq. (|10|l . 
where we have only one free phase factor ipi + ip2- Some 
further discussion is given in Appendix A. 

Note that the determination of the canonic state for 
the orbit specified by an arbitrary state vector of 
an n-qubit assembly can be formulated as a standard 
quantum control problem: the task is to find the global 
maximum of the vacuum state population given by the 



functional |(0| U\ ® . . . ® U n |\&)| 2 starting from the ini- 
tial state 1^). The space of the control parameters 
{P*, Pf , P*} is 3n-dimensional. Without taking advan- 
tage of additional structure, the complexity of this pro- 
cedure is in general exponential in n |61| . A possibility to 
improve the efficiency of this search is based on exploit- 
ing the solution of a set of differential equations which is 
discussed in Sect. IIII D 31 

We also note that, the requirement of maximum 
vacuum-state population alone is insufficient for deter- 
mination of the canonic state of assemblies consisting of 
qudits with d > 2. Indeed, a local unitary transformation 
not involving the vacuum state leaves this population in- 
tact, although it changes the amplitudes of other states. 
To eliminate such ambiguity, one needs to impose further 
constraints. As a possibility, one can maximize by a se- 
quence of step the populations of d— 1 states \k, . .., k), 
starting from k — and ending by k = d — 2. In this 
sequential procedure, maximization of the amplitude of 
the state \k, . . . , k) on the (k + l)-th step is done by a 
restricted class of local transformations belonging to the 
subgroup SU (d — k) that acts non-trivially only on the 
qudit states |m) with m > k. This algorithm leads to a 
generalization of the condition of Eq. ©: now, the am- 
plitudes of all states coupled to the states \k, . . . , k) by a 
single local transformation 6 SU(d — k) such as ?/>o...oi, 
V-'o...02, • ■ "01. ..12, • ■ ■) i>k...k,m>k but not xjJk...k,m<k van- 
ish. The action of these local transformations is indicated 
by dashed arrows in Fig. Hb). The remaining non-zero 
amplitudes normalized to unit vacuum amplitude char- 
acterize entanglement in the assembly of qudits fairly un- 
ambiguously. One has still to fix n{d— I) phase factors of 
unitary transformations, but in analogy with the qubit 
case, this can be done by setting real and positive some 
n(d— I) of d n — nd(d— 1)/2 — I non-vanishing amplitudes. 
In Sect. IV we discuss this in more detail. 



B. Nilpotent polynomials for entanglement 
characterization 

The amplitudes ip or a are not the most convenient 
quantities for characterizing entanglement, since they do 
not give an immediate idea about the entanglement struc- 
ture. For instance, for two unentangled qubit pairs, each 
of which is in a Bell state, one finds 

|*> = 5(|0,0) + |1,1))®(|0 1 0) + |1,1)) 

= ^oooo |0, 0,0,0) + ^ UO o 1 1,1, 0,0) 

+ Vooii |0, 0, 1,1) + Vim |1, 1,1,1), (11) 

where 

V'oooo = Vhioo = Voon = V-'im = 1/2 , (12) 
Quioo = c^oon = ami = 1 • 

In other words, though the system consists of two un- 
entangled parts, each of which is characterized by only 
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j> Since (cr^) — 0, these operators are nilpotent and they 

can be considered as nilpotent variables. Any quantum 
X\fO\\~ rea l \\flO\~ rea l \|/no"'" eQ / state |W) as in Eq. Q may be written in the form 



y i00=Q non> »-P,o,o> m=°- 



J22-— 202 



000 



I*) = O00...0 + ^oo...io-i" + • • • + ^01.. .oo-^-i + ^10. ..0^ 



| 100 > * |oio> * |ooi> +^oo...ii4-i + --- + ^ii...o ( T> r r_ 1 + ...)|0> 

»^|00t» = 2 ^ fe „-,.. fci nK + ) fc, |0>, (13) 

{fc;}=0,l i=l 

^) "220" ^ where each nilpotent monomial n"=i ( cr i + )' £t creates the 

■rjjr- 212 basis state \k n , . . . , k±) out of the vacuum state |0). Let 

210 \ ni 



F({ a i }) °e the nilpotent polynomial 



-m^4oo 022 ^L. 102 tw \ 

W V- 012^ = E fe^ n( ,^ (14) 



»W'~- * 002 {fei}=0,l 



containing only the zeroth and first powers of each vari- 
ed Piin= l able a t ■ A generic state |W) normalized to unit vacuum 

amplitude ipo...o = 1 can thus be written as F({cr^~})|0), 
B 1110-0 B 1101- B ion- Boiu- with F({af }) = 1 + a 00 ...ia+ + .. 

a fmo u moi> a |ion> H |om> Next; define the nilpotential /({cr+}) given by the log- 

R „,„,=! R,„,„™'»P'« R „,,„-! R,„m«»»P fa arlthm ° f F ({ CT J /)' 



BaaqSr,. fei 101> fw;^ r±<m<4 11{> Bjooi 

Hi°2°4ooo> B -=:foioo> &"£i°-&io> P-SSiljoooi 



/(K+})=ln[F({a+})] (15) 



P 0000 - l 



5ooo> {fci}=0,l i=l 



FIG. 1: a) Canonic form of the entangled state for three 
qubits. By a local transformation the amplitudes of the lowest 
excited states and the phases of the second highest exited 
states are set to zero. The amplitude of the lowest state is 
taken as the common factor determining the normalization 
and the global phase, b) In the case of qudits, the group 
of local transformations is richer, and the canonic form can 
be chosen such that it corresponds to maximum population 
of the symmetric states \k, fc, fc) where k = 0, . . . , d — 1. c) 
The structure of the tanglemeter for a 4-qubit system in the 
case of 51/(2, C) local transformations. The scaling factors qi 
can be chosen in many different ways, in particular such that 
some of the nonzero coefficients equal one. 

one parameter, three non-zero amplitudes a are present 
in the state vector. This is not convenient and a better 
description of the entanglement is desirable. 

We now introduce a technique which serves this pur- 
pose. Consider a standard raising operator 



acting in the 2-dimensional Hilbert space of the i-th 
qubit. Operators acting on different qubits commute. 



The coefficients /3fc Il / c „_i...fci and Oik n k n -\—ki can be ex- 
plicitly related to each other by expanding In F in a Tay- 
lor series around 1. This calculation requires at most 
n operations consisting of multiplications of the polyno- 
mial F — 1, which may generate an exponentially large 
(~ 2") number of terms. Note that /?oo...o = since 
c*oo...o = 1, so the nilpotential / starts with the first- 
order terms. Both F and / contain a finite number of 
nilpotent terms, at most 2" — 1, with the maximum-order 
term proportional to the monomial Yl7=i a t gi ven by the 
product of all the nilpotent variables. The canonic form 
of the state vector corresponds to a polynomial F Cl 

F c = l + a lJ <j+a+ ... , (16) 

which contains no linear monomials. The corresponding 
tanglemeter f c , also contains no linear terms and reads 

/c(K+}) = Pijafaj +..., (17) 

with Pij = a>ij. 

The discussion in Sect. Ill Al about the canonic form 
of the state vector applies to the tanglemeter as well. 
For most purposes, it suffices to employ the form which 
is unique up to phase changes of the nilpotent variables 
given by the local transformations er^ i— > cr+e 21 ^'' . The 
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coefficients ft therefore remain invariant up to n phase 
factors, unless these factors are specified by additional 
requirements. The phases may be chosen such that n 
of non-zero coefficients ft are set real and positive. Should 
the tanglcmctcr of the generic state be defined unambigu- 
ously, we can require this for the coefficients Pk n ...k 1 with 
X]i fcj = n — 1, in the same way as it was done for F. 
In special cases, where one or several such coefficients 
equal zero, some other conditions on the phases may be 
imposed. 

The tanglemeter / c ({Cj + }) immediately allows one to 
check whether two groups A and B of qubits are entan- 
gled or not. The following criterion holds: 

The entanglement criterion: The parts A and B 
of a binary partition of an assembly of n qubits are un- 
entangled iff 



d 2 fc({Xj}) 

dx k dx m 



, Vfc € A, Vm e B 



(18) 



In other words, the subsystems A and B of the partition 
are disentangled iff f A uB = /a({^6a}) + /slices}), and 
no cross terms are present in the tanglemeter. Note that 
the criterion of Eq. I|18|) holds not only for the tangleme- 
ter f c , but for the nilpotential / as well. However, we 
formulate the criterion in terms of f c , because the co- 



efficients P of the tanglemeter are uniquely defined by 
construction. 



C. Examples: Canonic forms for two, three, and 
four qubits 

For two qubits the result is immediate 

fc = Pw^t <r£, F c = 1 + a. xl a^a\ = 1 + /3ncr^ , 

(19) 

where the constant an = ft\\ can be chosen real. For 
three qubits, the canonic forms of F and / also differ 
only by the unity term, 

fc = (3 3 a+a++l3 5f 7+cT++f3 6 a+a++f3 7 cr+a+a+ =F c -l. 

.(20) 

Here, we have introduced a shorter notation by consider- 
ing the indices of ft as binary representation of decimal 
numbers, Oil — * 3, etc. One can make use of the fact 
that the variables o~f are defined up to phase factors, 
and set /?3, /?5, and Pq real. Expressing the invariants of 
Eqs. I©-© via the parameters in the canonic form, we 
obtain 



h = 1 + 2\p 7 \ 2 (/3f + (31 + Pi) + |/3 7 | 4 + 2Pl + Pi + 2plp 2 6 , 

h = i + m\ 2 (Pi + Pi + Pi) + \p7\ 4 + wi + Pi + wiPl , 

I 3 = 1 + 2\p 7 \ 2 (Pi + Pi + Pi) + \p 7 \ 4 + 2P 2 6 + Pi + 2PlPl , 
h + il 5 = 2 (Pi + 4p 3 p 5 p 6 ) . (21) 

This explicitly illustrates their linear independence. 
The tanglemeter for four qubits reads 

fc = Pso-to-t + P^t^f + pwi°i + pG°iai + Pw°to-+ + puo-t^ (22) 

+ P-io^o^oX + Pi3(ri(r£o-f + Pno-t a 2 a i + P\^t a t a 2 + P\^X a t a t a \ ^ 



while the coefficients cti of the polynomial F c differ from 
Pi only at the last position 

«15 = Pl5 + As/3l2 + P 5 PlO + P9P6- (23) 

One may note that the sums of the indices of the fac- 
tors in this expression are equal. The latter is a general 
feature for the relationship among the coefficients a and 
/?: the coefficients a are given by sums of terms, each 
of which contains a product of the coefficients P where 
the sum of the indices equals the index of a. We also 
note that a proper set of invariants of ®J =1 SUi{2) C* 
expressed via the components ip of the state vector can, 
in principle, be related to the tanglemeter coefficients, in 
analogy to the relation between Eqs. lloKil) and Eq. I|21|) 



for ®l =l SUi(2) invariants. 



D. Tanglemeter and entanglement classes for 
SL(2, C) operations 

We now take a larger class of local operations and con- 
sider arbitrary invertiblc linear transformations GL, in- 
stead of just the unitary transformations SU. Invertible 
transformations with non-zero determinant correspond in 
general to indirect measurements, that are measurements 
performed over an auxiliary system prepared in a certain 
fixed quantum state after it has interacted with the sys- 
tem under consideration. Besides allowing the realization 
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of measurements moregeneral than projective Von Neu- 
mann measurements Q, this procedure may serve as a 
tool for quantum control and quantum state engineering 
[62|. In the case where a single copy |(33| of a quantum 
state is considered, the outcome of the measurements is 
not achieved with certainty. Therefore, a stochastic fac- 
tor allowing for the outcome probability should be taken 
into account, whence the resulting state vector has to 
be renormalized in accordance. Since the normalization 
factor in the latter require an information about the ini- 
tial state vector, they do not strictly speaking form a 
group. In our approach, we do not impose a normaliza- 
tion condition on the wave function whatsoever and will 
be interested in finding the invariants [^3] of the trans- 
formations belonging to the group Q — <S>iSLi(2, C)(g>C*, 
where C* describes as before multiplication by an arbi- 
trary nonzero complex number and the transformations 
€ SXj(2, C) multiply the ith qubit state vector by a 2 x 2 
matrix of unit determinant. Another way to represent 
Q is to express it as the product ®iGLi{2, C) and fac- 
torize it over n — 1 redundant factors C*. The factors 
C* in each GLj(2,C) describe the same wave function 
transformations. 

We emphasize that the change of the wave func- 
tion renormalization does not result exclusively from the 
transformations S C*, but from some SL(2,C) trans- 
formations as well, and in particular, the transforma- 
tions ~ expqcrf with complex q. Thus, considering just 
(g>iSLi(2, C) instead of the full group Q may not have an 
explicit physical sense. Still, we will do it sometimes to 
better reveal the mathematical structure of the results 
obtained. Since the local SL transformations comprise 
the key part of Q 1 and following the established usus we 
will mainly refer to them and talk about sZ-entanglement 
in assemblies subject to indirect local measurements. 

Though local, SL operations can modify the set of 
quantities which characterizes entanglement in assem- 
blies subject to local SU transformations. Since SU C 
SL, many different su-orbits become equivalent under 
local SL transformations. In other words, the orbits of 
local SL transformations contain the su-orbits as sub- 
sets. Classification of si-orbits reveals the entanglement 
which persists despite the indirect measurements. In or- 
der to distinguish this type of entanglement from the 
invariants under local unitary transformations, one may 
call it sZ-entanglement. In particular, classification for 
three-qubit assembly |23| shows that all generic quan- 
tum states belong to one orbit of Q, which includes 
the canonic Greenberger-Horne-Zeilinger (GHZ) state 
(| 000) + 1 111)) A/2 The invariant Eq. © also known 
as 3-tangle r |2g , is different from zero only for this gen- 
eral orbit, and the value of r calculated for the state 
amplitudes normalized to unit probability descriminates 
different sit-orbits within this single general sZ-orbit. The 
states with r = 1 can be reduced to GHZ state by lo- 
cal unitary transformations, while for other states, with 
t < 1 indirect local measurements are required for the 
purpose. Moreover, there are five singular orbits of Q 



with r = that contain the states irreducible to GHZ- 
state. For four-qubit assemblies, the classification [24[ 
becomes much more involved, but still it gives an idea 
about the types of entanglement and eventual measures. 

Each element of the 5Lj(2,C) group, that is isomor- 
phic to the Lorentz group SO (3, 1), involves 6 parame- 
ters. In the general case, the number of invariants, 

D s i = 2™ +1 - 6n - 2 (24) 

is less than that for unitary transformations Eq. (121 . This 
counting is valid and returns a positive value for n > 4 
when the actions of different local operations are linearly 
independent. For n = 2, 3 where the number of the 
parameters in the group is more than the number of the 
parameters in the wave function, and the result of some 
local SL are redundant, no invariants exist. In particular, 
for two qubits, any generic state is equivalent under Q 
to the Bell state and for 3 qubits — to the GHZ state. 
For four qubits, there are 6 real invariants, for n = 5, 
D s i = 32, etc. 

A smaller number of sZ-invariants Eq. i|24|) as compared 
to that of su-invariants Eq. @ implies that different su- 
orbits may belong to the same sZ-orbit. In analogy to 
the su-canonic state, one has to define a sZ-canonic state 
as the marker of a sZ-orbit. In contrast to the unitary 
case where the canonic state has been defined by the 
condition of maximum reference state population, for SL 
transformations we introduce directly canonic form of the 
tanglemeter. To this end, we impose the following condi- 
tions: in addition to the requirement of the Eq. @, i.e. 
all n linear in a + terms of the nilpotential equal to zero, 
we require that all n terms of (n — l)-th order vanish as 
well. In other words, the sZ-tanglemeter takes the form 

n 

/c(K + })= £ ^...fallK^'. 

£jfci/{l,n-l} i=l 

and in this way we have specified 4n out of the 6n real pa- 
rameters of the local transformations that bring a given 
state to the sZ-canonic form. We thus left with 2n pa- 
rameters that have to be specified. 

In contrast to the unitary case, where the nilpotent 
variables a + are defined up to arbitrary phase factors, 
for SL transformations the variables in Eq. (|25|l are de- 
fined up to a complex-valued scaling factor af i— > 4li- 
One can further specify the sZ-tanglemeter by choosing 
these factors such that n complex coefficients of the tan- 
glemeter are set to unity. If convenient, one can impose 
another set of n requirements. 

As a first example, consider the three-qubit case. The 
sZ-tanglemeter Eq. (|25|l for a generic three-qubit state 
reads 

fc = 444, (26) 

where the coefficient is set to unity by the scale free- 
dom in the definition of the nilpotent variables. The 
corresponding wave function Fc is nothing but the GHZ 
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state. This shows again that all generic states belong to 
the same s^-orbit, which includes this state. There are, 
however, also three distinct singular classes of entangled 
states of measure zero whose tanglemeters do not in- 
volve the product o't' J 2 <7 3 anc ^ have one of the following 
forms, 

fc = <rt a t +< T 3 <T i > 

fc = O-t^X +< T t (T 2 > 

fc = <rt a x + a t a 2 ■ 

In this classification, we have only taken into account the 
states whose tanglemeters involve all three af such that 
no qubit is completely disentangled from the others. 

For a generic four-qubit state one finds the sl- 
tanglemeter 

fc = foet°Y + fc4°x + faUt + forfrt (27) 
+ (3 w al(j^ + {3\2ot a t + Px5^t a 3 a 2 a x > 

where the scaling factors of the variables af~ can be 
specified such that this form becomes equivalent to the 
expression given in Theorem 2 of Ref. |2J] : 



fc = As {<jfu+ + a+a+) + (a+a+ + a+a+) (28) 
+ & (a+a+ + <7+a+) 



cr, cr crq a A 



In Fig.njc), we illustrate the structure of s^-tanglemeter 
for this case with an alternative choice of the scaling fac- 
tors. 

It is worth mentioning that, though any generic nilpo- 
tential can be reduced to the canonic form of Eq. I|25() 
this turns out to be impossible for some sets of states of 
measure zero, as it is already the case for three qubits. 
These sets may play an important role for applications 
and can be grouped into special classes. Some of these 
classes are shown in Sect. IIII D 31 at the example of four 
qubits. There we also present an explicit algorithm for 
evaluation of s^-tanglemeters based on the stationary so- 
lutions of dynamic equations with feedbacks imposed on 
the parameters of local transformations. This yields the 
special entanglement classes in a natural way as singular 
stationary solutions. 

We conclude this section by discussing the precise 
mathematical meaning of the canonic states. The renor- 
malization of the wave function that follows the maxi- 
mization of the reference state amplitudes by local SU 
transformations belongs to the group C* of multiplica- 
tion by a complex number k. Therefore, strictly speak- 
ing, the applied transformations belong to the group 
®iSUi(2) <£> C*. However, the group ®iSUi(2) does 
not affect the normalization of the state vector, while 
the requirement ipo = 1 imposed on the canonic state 
uniquely specifies the number k thus allowing to intro- 
duce the tanglcmctcr as a characteristic of s^-orbits. In 
other words, once the condition ipo = 1 is satisfied, the 
group ®jS'{7i(2) ® C* becomes isomorphic to the group 
®iSJ7i(2). 



This is no longer the case for indirect measurements. 
Neither the group G nor its nontrivial part <S)jiSLj(2, C) 
conserve the state normalization. By imposing the re- 
quirement ipo = 1, we mark an orbit of G/C*, and 
thereby specify the structure of the canonic state given 
by the state amplitude ratios tpi/ipo expressed in terms 
of the s^-tanglemeter coefficients. However, a state of 
same structure but with a different normalization can be 
physically achieved in many different ways, - as a result 
of a single indirect measurement, or a sequence of two 
or more indirect measurements. The probability to ob- 
tain an outcome of the measurements that correspond 
to required G/C* transformation thus depends on the 
particular choice of the measurement procedure. There- 
fore, the complex factor n can be an arbitrary number, 
irrelevent to the values of the sZ-tanglemeter coefficients. 

However, when we consider just the nontrivial part 
®iSLi(2, C) of G, the factor n can bear certain physical 
significance. In fact, a transformation from ®iSLi(2, C) 
may bring a state initially normalized to unit probability 
to another one, which differs from the canonic state only 
by a factor k. In this case the factor k is uniquely defined 
function of the initial state 651. When the transforma- 



tion is unitary k amounts to 1/ yJ2i IV'il where the am- 
plitudes ipi of the canonic state are normalized to unity 
reference state amplitude, as required. For non-unitary 
SL transformations this quantity is different. There- 



fore, In 



can serve as a measure of non- 



unitarity of the transformation that discriminates differ- 
ent su-orbits that belong to the same G-orbit. 



E. How do the nilpotent polynomials relate to 
existing entanglement measures 



In general, there is no universal and precise definition 
of proper measures of entanglement 41], with the ex- 
ception of bipartite entanglement: as long as we are in- 
terested in entanglement between two parts A and B of 
a quantum system in a pure state, natural measures of 
such entanglement do exist. They are based on the re- 
duced density operator pa of either part, obtained by 
tracing over the quantum numbers corresponding to the 
other part B. In particular, S v n = — Tr[p^ log pa] and 
Si = 1 — Tr[p A ], give the von Neumann and the linear 
entropies, respectively 0], as already mentioned in the In- 
troduction. Clearly, both characteristics can be directly 
related to the tanglemeter parameters. However, the ex- 
plicit formulae giving these relations, which are simple 
for the case of two qubits 
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become awkward for larger numbers of qubits within the 
bipartition, as well as for higher-dimensional elements. 
This reflects the fact that the coefficients of nilpotent 
polynomials carry much more information about entan- 
glement than the simple bipartite correlations captured 
by the entropy measures. 

Another useful entanglement measure, concurrence C, 
has been introduced in Ref. 26] in the context of mixed 
two-qubit states, and has been employed for constructing 
the residual entanglement r, as a measure characterizing 
thr ee-q ubit pure-state entanglement and possibly beyond 
[2!i|6(j- Both C and r may be expressed in terms of the 
amplitudes tp of the sw-canonic state and in terms of 
the tanglcmctcr coefficients (3 Eq. (|20|l . The concurrence 
between the first and the second qubits reads 



2 |0ooo0no| 
2II&I- 



- |-0ioi^oii 

l&Asll 



■PI 



■Pi 



■Pi 



\p7 



(29) 



The residual entanglement, or 3-tangle, has the form of 
a fourth-order polynomial in the amplitudes. For the 
canonic state, it reads 



T3 



(0ooo0m) + 4 0ooo0no0ioi0on 
4|/3 r 2 +4/3 6 /3 5 /3 3 | 



(1 



Pi + Pi + Pi + \P 2 \Y 



(30) 



which up to a numerical factor is equal to the invariant 
\Ia + ih \ of Eqs. lf6l l21|l divided by the normalization i&c- 
tor£|V>| 2 = 1 + /3|+/3|+/3| + |/3 2 |. The presence of the 
normalization factor in the denominators of Eqs. (|29H . 

is due to the fact that these quantities are usually 
calculated for the state vector normalized to unity while 
the coefficients (3 refer to the tanglemeter that is the log- 
arithm of the canonical wave function with the normal- 
ization "0000 — 1- 

What are convenient measures that can be introduced 
to characterize s^-entanglement? We have seen that all 
generic states of the assembly of three qubits belong 
to the same orbit of Q and strictly speaking there are 



no invariant measures at all. However, the sw-invariant 
I4 + H5 of Eq. @ remains invariant under the restricted 
class of transformations ^ =1 SLi(2, C), while the other 
SU invariants ii 2,3 of Eq. (JSJ) depending on both 
and 0* change under SL transformations. Hence, in 
this restricted sense it may serve as a measure for sl- 
entanglement. 

The measures characterizing the s^-entanglement for a 
4-qubit assembly can be constructed in a similar way. We 
take products of several factors ~ (but not the factors 
~ tp*) and convolute it over SU(2) indices with invariant 
tensors e M |66|. The simplest combination 



i (2) = Ajkii> m 



(31) 



is si-invariant and can be taken as a characteristic of sl- 
cntanglement, remaining not invariant only with respect 
to the transformations S C*. There are three different 
sWnvariants ~ tp 4 , 



r(4) _ r(4) 
J 12 — J 34 



(4) 
13 

(4) 
14 



-(4) 
24 

-(4) 
23 



tP lk3 l^ mJn tPo mp ni> 0kpl , 



(32) 



are 



The ratios 7^/(/( 2 )) 2 , I^/{I^f, and I&7(I< 2 >) 
in addition invariant with respect to multiplication of 
the state vector by an arbitrary complex constant and 
thereby they are invariants of Q. Were these ratio lin- 
early independent, they would give us a complete char- 
acterization of 4-qubit entanglement, since the 4-qubit 
s^-tanglemeter Eq. (|28(l involves 3 complex parameters. 
However, they are not. The following identity 



r( 4 ) 



r(4) 
'13 



(33) 



makes these quantities inconvenient for the entanglement 
characterization. 

We therefore turn to the 6-th order invariants and con- 
sider following three functionally independent combina- 
tions 



I\2 q (jPingd'Pmrko^Psjph 'Pingo'Pmrkh'Psjpd) ^ tp ^ tp ^ , 

-^23 g (jPijpo'Pmngh'Psrkd 'Pijpd'Pmngo'^Psrkh) ^ ^ tp ^ , 

J 13 J = I WijkhTpmnpdlPsrgo ~ ^ijgh^mnkd^srpo) ^d^nph^sjko ^ (34) 



whose differences give the invariants Eq. (|32|l multiplied state is awkward. However they take a simple form 
by I^ 2 \ Explicit form of these invariants for a generic 

I (2) =0oooo (t + x + y + z) 
III = 4 0oooo (t + x-y-z)(tx- yz) , 
4t = 4 0oooo (* ~ x + y - z ) {ty ~ xz ) > 
l[f = 4 0o 6 ooo (* ~ x - y + z) (tz - xy) , (35) 
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for the canonic state, where the sZ-tanglemeter Eq. 128(1 
suggests 



V'llOO 


= ipoon = 


V'ooooV^ 


V'lOOl 


= "01001 = 


00000 \Fu 


V'OlOl 


= ^1010 = 


00000 \fz 


ipun 


= 00000^ 





for the nonvanishing amplitudes, and where the notations 
x = 0\, V = 0%, z = 0l, and t = ip nn = 1 + x + y + z 
are employed. We introduce new variables 

x = 0oooo (* + x ~ V - z ) , 

Y = 0oooo (t-x + y- z) , 

z = 0oooo (t-x-y + z) , (36) 

and find 

l[f = X{I^X-YZ), 
4f = Y(I^Y-XZ), 
l[f =Z(lWZ-XY). 



Solving this system of equations yields 

X = y/(l[*>+P)/lW , 

Z = yJ{lg + Pyi(*), (37) 
where P is a root of a cubic equation 

Different roots of these equations and different signs of 
the square roots in Eq. Q37JI yield different sZ-canonic 
states, that either coincide or are related by SL transfor- 
mations. The amplitudes of these states can be written 
explicitly 



0000 



ig + P+Jig + P + Jig + P- (mf 2 



y/2{lW) 

01100 = 00011 = - 



1/4 



I®+P-Jlg + P-Jl® + P + (l(™ 3 '> 



2(K 2 



»V4 



•01001 = "01001 



2(K 2 )) 



1/4 



00101 = V'lOlO = 



I^ + P-Jj^ + P-Jl^ + P + (I^ 



2(K 2 



,V4 



0im = V j V V (38) 



2V2 (K 2 >) 1/4 J Jlif +P+ Jl® +P + Jl$ + P- (K 2 )) 



,3/2 



while the ratios V'noo/^oooo: ''Aiooi/V'oooo, an d the 4-qubit assembly can be completely characterized by 
Vtoioi/^oooo yield the sZ-tanglemeter cocfhcicnts three independent scale-invariant complex ratios 

05, and 06, respectively. Thus, the sZ-entanglement in 
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03 



05 = 



06 



ig+P-y/lg+P+iX 



(6) 



,3/2 



V2y 


' /r(6) 
V J 13 


+ P 4 


1 r(6) 
- V 7 23 


+ P4- 


- V J 12 


+ p 


- ( /(2) ) 


\J\ 


+ 


P- 


A ( s 5 + 


P- < 


+ 


p+ 


(/(2)) 3/2 






+ P4 


- A ( 3 3 


+ P4- 


- ^? 


+p 


- (/(2)) 3/2 




/4? + 


P- 


+ 


P- < 


As 5 + 


p+ 


(/(2)) 3/2 



WW + p + W + ^ + + p - 



coming from the invariants Eqs. (|34I31I) . 

As for a measure characterizing entanglement for 4 
qubits, one has to consider at least two quantities. The 
first is the sum \ip 2 \ over the amplitudes Eq. (J3SJ, 
which gives the regular normalization of the canonic-like 
state. Once the invariants of Eqs. (pTTJl . ffity are calcu- 
lated for a state normalized to 1, this sum shows how 
the SL transformation required for setting the state to 
the canonic form is different from a unitary transforma- 
tion, whence | In ^ 1 4> 2 \ | provides us with a measure of 
this nonunitary. The root P and the signs of the square 
roots Eq. I|37|) have to be chosen such that |lnj] \ l l> 2 \\ 
is minimum. This quantity discriminates different su- 
orbits that belong to the same sZ-orbit in analogy to the 
3-tangle, which discriminates different su-orbits within a 
single generic SL orbit of 3-qubit assembly. It may serve 
as a measure of su-entanglement within a sZ-orbit. The 
second quantity has to discriminate different sZ-orbits 
and serve as a measure of sZ-entanglement. A natural 
candidate for that is the sum of moduli squared of the sl- 
tanglemeter coefficients j3, which takes value for GHZ 
canonic state, remaining larger for all other states. The 
choice of P and the signs in Eq. 1|37[) has to be done such 
that this quantity is minimum. This measure shows us 
how close is the orbit to the GHZ orbit. 

One may also ask for simpler measures that would 
be polynomials on the state amplitudes V- It 
turns out that two such characteristics, \l^\ and 



(6) 
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r(6) 



'23 



r(6) 



! 12 



) / l-^ 2 ^ 2 ! can be directly associ- 
ated with the su and sZ-measures. In Fig. [21 we show 
these characteristics plotted versus exp | In | -0 2 1 1 and 
X]|/3 2 |, respectively, for a variety of ~ 10 2 randomly 
chosen assembly states normalized to unity. One sees 
that |P 2 ^| strongly correlates with the measure of non- 
unitarity, while the sum of the moduli of 6-th order in- 
variants majorates the sZ-entanglement measure based 
on the sZ-tanglemeter coefficients. Other combinations 
of invariants do not correlate with the tanglemeter coef- 
ficients. 



III. QUANTUM STATE OPERATIONS AND 
DYNAMICS IN TERMS OF THE NILPOTENT 
POLYNOMIALS 

In this section, we first describe the effects of local and 
gate transformations on the assembly state vector as al- 
gebraic manipulations of the corresponding polynomials 
F and /. In principle, by applying a properly chosen se- 
quence of finite local transformations, one can reduce a 
nilpotent polynomial to the canonic form, thereby spec- 
ifying the tanglemeter. However, straightforwardly ap- 
plying these transformations is not a very practical way 
to proceed, since it usually requires lengthy calculations. 

We therefore turn to infinitesimal transformations sec- 
ond, and derive the equations of motion describing the 
dynamics of the nilpotential under continuous local and 
gate operations. We show that, for an important class 
of Hamiltonians supporting universal quantum computa- 
tion, the dynamic equation for the nilpotential acquires 
a well-known Hamilton- Jacobi form. 

We thirdly demonstrate how to determine the tan- 
glemeter with the help of such equation. To this end, 
a proper feedback is required, ensuring that the param- 
eters of infinitesimal SU(2) or SL(2, C) transformations 
are adjusted to track current values of the nilpotential co- 
efficients. The tanglemeter appears as a stable stationary 
solution, that is a focus of the resulting equation. We il- 
lustrate this method in the example of four qubits. We 
show how to find the sZ-tanglemeter for a generic 4-qubit 
state and how to explicitly identify a number of special 
classes that cannot be reduced to this form. For these 
classes we suggest alternative natural tanglemeters. 



A. Local operations 

A general local unitary transformation Eq. @ applied 
to i-th qubit can also be expressed in the equivalent form 

U i = e A *'*e Bia *e c " r *, (39) 
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a) 



0.6 



Oa 



I 



(2) 



The explicit expressions 



A, = 



(iPf-Pf)sinP 
PcosP + iP/sinP' 



B t = log 



cos P 



iP 2 sin P 



P 



= (iPf + Pf)sinP 
' P cos P + iP/ sin P ' 



£ Of?) 



(40) 



b) 



0J 



exp\ln 



1.2 



1.4 



/ • .» • 

/'• -»...:.. ••. 

_, , t , I ' , , I , L 

0.1 



i (6 Lt (6 L t (6) 



Q(2))2 



0.2 



relate the parameters in Eq. and Eq. 

The transformations Eq. (|39(l act on the state vec- 
tor \^f) — F({af})\0) and yield a transformed state 
F'({af})\0). One can formalize the rules allowing one 
to obtain F' from P. Bearing in mind that <t _ |0) = and 
cr z |0) = — 10), one can represent the action of af , of , a~ 
as appropriate differential operations for the nilpotent 
variable af. The application of the operator af is 
straightforward - it is a direct multiplication: this oper- 
ation eliminates the terms that were proportional to af 
prior to the multiplication. The application of o~ is a 
kind of inverse: it can be considered as a derivative with 
respect to the variable af, which eliminates the terms 
independent of af and makes the terms linear in af in- 
dependent of this variable [63] . Finally, the application 
of af changes the signs of the terms independent of af , 
and leaves intact terms linear in af . These actions are 
summarized by the following formulae 



FIG. 2: a) The polynomial invariant plotted versus the 

non-unitarity measure exp |ln^ |"i/> 2 | | for a set of ~ 10 2 ran- 
domly chosen assembly states, b) Same for the combination 

7 i2'|) / | /(2) | and the E |/3 2 |given by the 



Hi) 



si-tanglemeter. 



afF = a+F, 
OF 



vfF 



da: 



F + 2af 



dF 



1 dal 



(41) 



which better suits the consideration in terms of nilpoten- 
tial polynomials, since each step of the local transforma- 
tion can be expressed as an operation linear in er, that 
is 



1 



Ciof , 



cos (Bi) + sin (Bi) a\ 
1 + Ai<7- . 



while each unitary operation U (af , cr? , of) 
can be represented by a differential operator 



By sequentially applying the three transformations of 
Eq. Ij39(l to F, a local transformation can be interpreted 
as multiplication by an exponential function of af , fol- 



lowed by a linear transformation ay 



J2Bi 



")of 



the variable of and multiplication by e ' , leading to 



UiF(af,...,af, ...,al] 



~e B ^F{af,. 
G(af,...,e 



2 LI, 



e 2i V+ 



~e B ^G(af,. 



(42) 



Since local operations on different qubits commute, this generalized to n qubits. 

single-qubit transformation may be straightforwardly Note that in order to cast F into the canonic form, 
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{fc 4 }=0,l i=l 



{fei}=0,l;fci+...+fc n /l t=l 



one has to solve a set of nonlinear equations for the pa- 
rameters Ai, Bi, and Cj. This can be done explicitly 
only for at most four qubits, while for a larger system an 
efficient numerical technique is required. This task can 
be accomplished by an iterative procedure in the spirit of 
the Newton algorithm, that is, by consecutively applying 
a series U n . . . U-zUi of linear transformations {/,, each of 
which eliminates the terms linear in af ■ However, this 
procedure may require infinitely many iterations, since a 
linear transformation applied to one of the o~f may (and 
usually does) generate terms linear in other crTV . In 
Sect, im D 31 we show how dynamic equations describing 
the evolution of the nilpotential / under local transfor- 
mations offer a better tool to solve this problem. 



The nilpotent polynomials Ai — Ai ^cj^^- H and 

Aj = Aj (( j }) wn i cn are the coefficients in front of 

the variables of and of , respectively, undergo a unitary 
rotation 

Aicrf i— > Aicrf cos t + ^Ajaf sin t 

Ajcrj h-> Aj-crt cos t + L4;o-+ sin t . (44) 

in the same way as the components of a qubit state vector 
do under an SU{2) transformation. 



B. Two-qubit gate operations 

Quantum gates are unitary transformations acting on 
finite subsets of qubits in the assembly. In particu- 
lar, two-qubit gates Uij operate non-trivially on the pair 
Thanks to general universality results 0iE3j an 
arbitrary non-local transformation on n qubits may be 
expressed as a finite sequence of arbitrary single-qubit 
operations and two-qubits operations drawn from a stan- 
dard set, applied to both individual and pairs of qubits 
according to a certain quantum network. Thus, start- 
ing from an initial computational state, any state may 
be reached through the application of a quantum circuit 
built from gates in the set. We consider here the simplest 
choice for the standard two-qubit gate operation, 



Uij = exp [ii(a+<r- + o i rr+)] = exp 



9 ^ y ■> . 9 " 

= cos - + o,o~, sin - 

9 ' J 9 



sint, (43) 



depending on the single parameter t G M, where the ten- 
sor product symbol, of a* — of ® o* , is implicit. 

Only the terms of F that contain af or cr+ are affected 
by the transformation Eq. I|43|l . The terms that either do 
not contain these variables or are proportional to their 
product are left intact. 



C. Local and gate operations in terms of the 
nilpotential 

Equations (|42|I - H44|) are particular cases of general ex- 
pressions for transformations of nilpotent polynomials F 
under the action of unitary operations. We now consider 
this in terms of the nilpotential. 

The rules of Eq. (|39|) , allows one to express the action 
of a unitary operation U ({ erf , o~ , erf } ) as a differential 
operator acting on the nilpotent polynomial F 

p=u {{"''it M ^f- 1 }) F ' (45) 

while for the nilpotential one finds 



f'=log[u[\o. 



1 ' do+ 



d n + d 

2cr+ 



0, 



o 



(46) 



Note that a generic transformation Eq. H46fl of an ini- 
tially canonic polynomial does not necessarily results in 
another canonic polynomial. 

Let us consider two particular cases of the general 
transformation of Eq. (|46|l : (i) A local unitary opera- 
tion Eq. @ with Pf = Pcos0, Pf = Psin^, P z = 0; 
(ii) The two-qubit gate Eq. 143(1 . They transform the 
nilpotential according to 
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and 



f = f + 2afaf (1-cost) 



a 2 / 



[aj — aj cos t 



sin 2i 



• + • ^ 9f 
a a J 



af cos t — i af sin i) 



a/ V . . a/ a/ , sm2t ( df 

1 + 2i sm i 



da? 



df_ 

daf 

"\ (48) 



respectively. 



D. Equations of motion for the nilpotential 

Consider now an infinitesimal unitary transformation 
U = 1 — i dt H which is not necessarily local. The incre- 
ment A/ of the nilpotential / suggested by the Eq. 147|) 
reads 



A/ = log (Ue f ) - log (e / ) = log (l - i dt e~ f He f ) 

This yields the following dynamic equation for /, 

■ 9 f -fu f 
i— = e J He J , 

at 

which wc discuss in detail in the rest of this section. 



(49) 



(50) 



where Pf = Pf ± iPf, and we first separately consider 
only the term H l = Pr(t)af + Pf(t)a7 + P?{t)af in 
the sum. Upon substituting it in Eq. I|5U|) and splitting 
the nilpotential / on the right hand side in two parts, 
the part /o = f — afdf /daf independent of af, and the 
part /i = afdf/daf linear in af , we obtain 



= e -*t9f/d°t H^tdf/d.t . (52 ) 



1. Local Hamiltonians 



We begin with the case of a local Hamiltonian 
H = J2 + P?{t)af + P*(t)ai 

i 

= ]T Pr(t)af + Pf(t)ar + P*(t)a* , 



The part fo commutes with the derivatives entering 
the Hamiltonian and therefore cancels. Substitution of 
(51) Eq. (gTJ into the Hamiltonian of Eq. J^U followed by 
expansion over the nilpotent variable af results in 

I 



i^- = -P z + p-a^ 
dt 



Of 



^Pfaf+P, , 

oa 



df_ 

daf 



Straightforward generalization of this equation to the case of the Hamiltonian Eq. (|51|) yields 

()t =it\- p * + rr*t + (ip?4 + p+) — - p+ - + ( — 2 



.df ^ 

i=l 



daf 1 1 \da. 
Another equivalent form of the same equation reads 

{ % = E l p ? - 1) + Pf (fi + 4 fM) + (A 4 - ffr 

i—l 

i 



(53) 



(54) 



(55) 



where we denote fi = df /daf . Note that the coeffi- the right hand-side of Eqs. (|53[l - i|55|) does not depend on 
cients p*- y ' z can be functions of time. Also, note that 
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the constant term in /. In other words, the latter, though 
evolving with time by itself, does not affect the evolution 
of the "essential" coefficients in the nilpotential in front 
of the nilpotent variables and their products. 



2. Binary interactions 
We now consider the binary interaction 
H = E Gy^K ?' k,x=+,-,z, (56) 

among the qubits. Note that the local transformations 
Eq. (|55|l can be absorbed into the time dependence of 
the coupling coefficients G^{t) by simply passing to the 
interaction representation. In order to achieve universal 
evolution in this representation, one needs to consider 
all nine coefficients G^* characterizing the interaction of 
Eq. (|56[) between a pair {i,j} of qubits as being different 
from zero. An alternative way is to chose such a repre- 
sentation that tensor G^{t) in Eq. 1(56(1 takes the form 
of a diagonal spherical tensor with respect to the upper 
indices. In this representation, the Hamiltonian 

h = j2 p?m + E G r (*) fa + °7 + °r<tf) (57) 

i,K i,j 

+ ^Gf/(tKa}+^G++(t)(a+a-++<7rar) , 



of motion for /. It reads 

(59) 

%n 1 

Note that Eq. 1(59(1 formally resembles the Hamilton- 
Jacobi equation for the mechanical action of classical sys- 
tems with the Hamiltonian 

H = E [ p r(t>i + p t(t)Pi (i - xaH)] 

i 

+ Gij (*) [xjPi (1 " %iPi)] . (60) 

where pi = df /daf plays a role of the momentum, while 
Xi = (T 2 + are the coordinates. Comparing with the con- 
ventional classical Hamilton- Jacobi equation, the only 
essential difference is the factor i multiplying the time 
derivative and the presence of complex parameters that 
can be interpreted as time-dependent forces and masses. 
After cumbersome calculations taking into account the 
fact that the constants of motions entering the action 
function are nilpotent variables, one can reproduce the 
finite transformations of Eqs. 1(47148(1 . 

3. Dynamic equations for the nilpotential and construction 
0/5(7(2) and SL(2, C) tanglemeters 



apart of the local operations Eq. (|54|l involves also the bi- 
nary interactions determined by only three real coupling 
parameters Gff(t), G±~ (t) = G^+(t), and G+ + {t) = 

G--{t). 

The explicit forms of the equations of motion Eq. I|5U|) 
for the Hamiltonians Eqs. I|56 |) .l|57 |) are rather awkward. 
We note, however, that universal evolution is achieved 
|55| with an even simpler Hamiltonian 



H^P+^at+^Pr^ar (58) 

i i 



with -Pj + (i) = P~(t)*, which depends on a smaller set 
of operators, af, o~ , and (af <rj + cr~o-j~) . Repeated 
commutators of these operators satisfy the Lie-algebraic 
bracket generation condition for complete controllabil- 
ity, that is, all-order commutators span the full space of 
Hermitian operators for the assembly, and thus ensures 
universal evolution. It therefore suffices to specify the 
form of Eq. (JSUJ for the Hamiltonian of Eq. (JSSJ. 

In Appendix C, we derive the corresponding equation 



The dynamic equation Eq. 1)59(1 suggests an algorithm 
for evaluation of the tanglemeter. This is based on the 
idea of feedback by adjusting the parameters Pi of the 
local transformations in Eq. (|59|l in function of the cur- 
rent values of the tanglemeter's coefficients. To this end, 
we fix the terms linear in in the local Hamiltonian of 
Eq. JSU as 



pr 



where 



Pi 



Of 



da 



(61) 



(62) 



are the coefficients of the linear terms in the nilpotential 
at a given time. From Eqs. I|59I61|I we find the evolution 
of these coefficients under local transformations, 



d 2 f 



3 da+daj 



- ft IA 



(63) 



When / is close to f c , the matrix of the second deriva- 
tives Mij = (<9 2 / / 'd&j dal ~)\ a >o can be explicitly ex- 
pressed in terms of the "second excited state amplitudes" , 
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^o...fei=i,fe,=i,...0) which enter Eq. (JSJ and were introduced 
when discussing the su-canonic form of states. According 
to Eq. (EJ, 



11 da+d*+ 



=Pij = CXij = ^O...k < =l,*j=l,...o/V'0,...0 • 

The condition of the maximum reference state population 
for a state in the canonic form suggested by Eqs. i|8l9|) 
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FIG. 3: Diagram illustrating the dynamics of the coeffi- 
cients j3i of the nilpotential / for a 4-qubit system subject 
to 51/(2, C) transformations. The coefficients corresponding 
to quadric, cubic, and bilinear terms are shown by black, 
gray, and white circles, respectively, a) Flux directions. The 
SL(2, C) transformations are applied in such a way, that the 
polynomials remain in the canonic form with respect to su(2) 
transformations, with all linear terms Pi<jf = 0. The dashed 
arrows show the linear contributions of different states to the 
time derivatives of the neighboring states, and the dotted ar- 
rows jointed by asterisk depict two types of bilinear contribu- 
tions. The other, similar terms may be constructed by sym- 
metry, b) Canonic form for the generic state and for singular 
classes corresponding to one c), two d), three e), and four f) 
vanishing eigenvalues of the determinant Eq. I|ri9|l . 





2 i 




- 1 



implies that the population increment, 

= -Re | Vo...o| 2 X 9t (f^+^Mij + 9j , 

is always negative. As the phases <fi are arbitrary, it 
follows that the eigenvalues of My lie within the unit 
circle and hence their real parts lie in the interval (—1, 1). 
Therefore the equation (|63f) linearized in the vicinity of 
the canonic state 



m 

dt 

m 



(64) 



has a stable stationary point at {Pi} — 0, which implies 
that the coefficients j3i standing in front of the linear 
terms (3i(jf in the nilpotential tend to zero exponentially. 
The presence of the nonlinear terms (3i |/?i| 2 yet accel- 
erates this trend. Therefore, an arbitrary nilpotential 
/ subject to the local transformation with the param- 
eters of the Hamiltonian Eq. (|51(1 chosen according to 
the feedback conditions Eqs. Hbllb2|) . rapidly converges 
to the tanglemeter f c . The problem of finding an effi- 
cient numerical algorithm for determining the tangleme- 
ter for large assemblies is thereby solved. Verification 
that the outcome indeed corresponds to the global maxi- 
mum of the reference state population should finalize the 
procedure. Note however, that the maximum vacuum 
state population obtained with local transformations cor- 
responds to the maximum population of the ground state 
for each qubit. On the other hand, for a given set of 
single-qubit density matrices, the local operations max- 
imizing the ground state population of each qubits are 
uniquely defined. Therefore, the only maximum of the 
reference state population is the global one, and hence, 
no matter what the initial state is, the procedure indeed 
converges to the canonic state and no verification is re- 
quired. 

A procedure of reducing the nilpotential to the canonic 
form can be carried out also for SL transformations. 
At the first stage of this procedure, we reduce it to the 
su-canonic form so that the terms linear in o1[ are ab- 
sent. Then we apply SL operations. An element of the 
5X(2,C) group can be represented as exp{it(P~cr i f + 
P^a^ + P?af)}, where Pf and are no longer com- 
plex conjugates and Pf is also a complex number. 

Finding the sZ-canonic state can also be formulated as 
a control problem, based on the feedback. We choose 
the parameters Pi in the Hamiltonian Eq. I|54|) in such 
a way that the terms in the nilpotential involving the 
monomials of order one and of order n — 1 in <jf would 
decrease exponentially with time. To this end, we may 
choose at this stage P[ = and impose two conditions, 
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(i) the condition 

n 



i=l 



d 2 f 
da+da+ 



expressing P- via P^~ which is keeping the nilpotential 
in the form of tanglemeter, and (ii) the condition 



four third-order coefficients to zero, we gradually reduce 
the SU(2) canonic form of Eq. (23 to the SX(2,C) form 
of Eq. (|2"7|) . This control process results in an exponen- 

= -l^Khj, (65) 

tially fast vanishing of the coefficients 07, 0i3, 0n, and 

/3l4- 

We now write down explicitly the differential equations 
for these coefficients, 
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i/5i4 




- 2P+0 6 W + 2P+0 6 12 


+ 2P+0 lo 0i 2 




= 2P+/? 5 /? 9 


- P+p 15 + 2P+j3 5 p n + 


2P^ 09012, 




i/3n 


- 2P+/3 3 /3 9 


+ 2P+j3 3 w - P+0 15 + 


2P 4 09010, 




i/3 7 


= 2P+/3 3 /3 5 


+ 2P+0 3 6 + 2P+0 5 6 


- P 4 + 015 , 



which ensures the exponential decrease of all n coeffi- 
cients in front of the second-highest order terms. 

After having eliminated the monomials of orders f and 
n — l, we can specify the scaling parameters Pf such that 
n additional conditions are imposed on the tangemeter 
coefficients. For example, one can set to unity the coeffi- 
cients in front of the highest order term and set (n—l) 
coefficients in front of certain monomials equal to (n—l) 
coefficients of other monomials. Within the group G, 
multiplication by a complex number k G C* allows one 
to normalize the canonic state to unit vacuum-state am- 
plitude and thereby to get rid of the constant term in the 
sZ-tanglemeter fc Eq. lj2~5l) . 

The condition Eq. I|66|l on Pl~ is written implicitly as a 
set of n linear equations. These equations can be resolved 
for generic states as we will show in the next section in 
the four-qubit example Eq. i|67[l . However, they have no 
solution when the determinant of the system vanishes. 
These singularities correspond to singular classes of en- 
tangled states and require special consideration. 



4- Example: Classes and si -tanglemeter s for 4 qubits 

Now with the examble of 4-qubit assembly we illus- 
trate the procedure of evaluation of the sZ-tanglemeter 
with the help of dynamic equations supplimcnted by the 
feedback conditions. The sw-tanglemeter Eq. (1221) has 
11 complex coefficients 0i and one has to solve a sys- 
tem of eleven first-order nonlinear differential equations. 
Instead of presenting this awkward system explicitly, in 
Fig. |3| we schematically depict contributions to the time 
derivatives of the 0i which are either linear or bilinear 
in the tanglemeter coefficients, and can be interpreted 
as a sort of "entanglement fluxes" [2fT |. One notices 
that the coupling of the second-order terms 0ijaf<r^ 

to the fourth-order term 0i5af a 3 f a 2 <Ji occur via the 
third-order terms 07<r 3 a 2 a^, P13&4 o^ cr^ 1 0ii a 4 °2 °2 > 
0i4(J 2 a 3 <j\ , and thus the time evolution of all 0i stops 
when this third-order coefficients 07, 0i3, 0, and 0u van- 
ish. Therefore, by setting the time dependence of the 
parameters P^~ , P 2 , P^ an d P4 such that they drive all 



(67) 



and see that, in the general case, by a proper choice of the 
paremeters P^~ one can impose the feedback conditions 
such that these equations take the form 



7 = -0T, 0u = -0n; 



13 



-013; 



14 



-01 



The evolution implied by these equations brings the 
nilpotential in the form of sZ-tanglemeter fc- We note 
that the coefficients P~ should satisfy the requirement 
Eq. (|65|l which ensures that the nilpotential always re- 
mains in the form of the tanglemeter f c during this evo- 
lution even if the state does not remain in the same su- 
orbit. We thus arrive at the tanglemeter fc of Eq. 127fl . 
defined up to the scaling factors. Now we can invoke the 
scaling of the nilpotent variables errand reduce the tan- 
glemeter f c to the sZ-canonic form fc of Eq. Q28j l , unless 
one of the bilinear coefficients vanish. The latter case 
corresponds to a measure-zero manifold and the canonic 
form may be chosen as 



fc = <rt<4 + 05 (4°! + 



2 , (68) 

06 {vtcri + °t °2~) + a 4 a t a 2 °1~ > 



or in any equivalent form resulting from a the permuta- 
tion of the indices. More singular classes are discussed in 
Appendix B. 

Reducing / to the canonic forms of Eqs. I|28I68I) is 
unattainable when the determinant 



-01 5 20 6 1O 20 6 0i 2 20^012 

20509 -015 2050U 20 9 12 

20 3 9 20 3 1O -0i5 20 9 1O 

20305 2030 6 20506 -015 



(69) 



vanishes, and it becomes impossible to impose the re- 
quired feedback conditions. In that event, we loose con- 
trol over the dynamics of 07, 0i3, 0ii„ and 0i4, and some 
linear combinations of these coefficients cannot be set to 
zero by a proper choice of P^, Consider this singular 
case in more detail. The determinant Eq. I|69|l is equal 
to zero when one or more of its eigenvalues, 
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71 = 0i 5 - 2^050609010 + 2^030609012 - 2^0305010012 , 

72 = Pis + 2^050609010 - 2^030609012 - 2v/0 3 050io0i2 , 

73 = /?i5 - 2v/0 5 6 9 lo - 2^/030609012 + 2 ^0305010012 , 

74 = 0i5 + 2 v //? 5 A i /39/3io + 2^030609012 + 2^0305010012 , (70) 



vanish. Let us first focus on the case where only the first 
eigenvalue is zero. This implies that six coefficients in 
front of the bilinear terms and the coefficient 0i5 in front 
of the 4-order term are no longer independent parameters 
- the last one being the function of the first ones explicitly 
given by 7i =0. The eigenvector 

(- V0601O012, V0509012, - V030901O, V foPsfo) > 

corresponding to 71 gives the combination of the cubic 
I 



terms that cannot be eliminated, 



- W 0601O012O- 4 f <7.+ (7+ + A\/0509012CT4 <7+0- 4 



Clearly, this combination is determined up to a scaling 
factor A. The nilpotential of Eq. ill' 2 1) thus takes the form 



/ = 3 cr^cr ] f + 05^3 a X + @9 a 't a l + G t + 010 04" °2~ + Pl2 a t CT 3 

+ A ^030506^0-^0-+ - ^03090100-^ 0-+0-+ + v /0 5 090i2Cr 4 h CT + cr + - yj 0601O012O-4 &t a 2 ) 

+ 2 (V05060901O ~ ^030609012 + V030501O01 2 ) 0"+ a+<7+a+ , (71) 



specified in terms of seven complex parameters (A,0i). 
We can eliminate four of these complex parameters by 
setting the scaling factors of (7+ , and arrive at the form 

f C = 03 + vUt) + 05 (4^1 + + °M) 

+ 6 ((7+0-+ + 0+0+) 

+ 0+0+0+ — 0+0+0+ + 0+0+0+ — 0+0+ 0+ 

+ 2 (0 5 6 - 030 6 + 030 5 ) CT+O+O+O+ > ( 72 ) 

which depends only on three complex parameters. This 
combination can be considered as a class of the poly- 
nomials that cannot be reduced to the canonic forms of 
Eqs. (|28I68[1 by the sequential application of infinitesi- 



mal transformations preserving the canonic SU (2) form. 
Permutation of indices of of give equivalent classes. 

Next, we consider the case where two of the eigenval- 
ues, say 71 and 72, of Eq. (jTOjl are zero, that is, 



015 = 203012 , 05010 = 03012 ■ 

The corresponding eigenvectors 

("V0601O012, ^0509012,-^0309010, ^030506 
("V0601O012, -V0509012, V030901O, ^030506 

suggest the form of the nilpotential 



J 



/ = /3 3 o+o+ + /3 5 o+ 0+ + (i 9 o+o+ + /3 6 o+o+ + P w o+o+ + 0i 2 cr 4 ;<7 4 " 

+ A (V030506O-3V+O-+ - \/030901Off4^2 f O- + + ^05090120^ O^ 0+ - ^06010012^ o£o+^J 
+ M (\/030506 0-^(7+ (7+ + ^03 09010(7^ <7 + (7+ - V05090i2CT 4 f O-3 f O- + - ^0 6 01 01 2 44^) 
+ 203012(7+ f7+CT+ ( 7+, 



which after a proper scaling can be simplified Again, this depends on three complex parameters and 

fC = 0't <7 l + ^i 17 ! + °X a 2 + °t a t + 06 (04" (7+ + 0-+CT+ ) 

+ 07 (0- 3 f 0-2 f ai f - o-|cr+(r^) + 0ii (0-^02" cr + - ctI ^+(7+) 
+ 2cr+ 0-+0-+ <7+ . (73) 
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leads to equivalent classes under indices permutations. 
Similarly, the case where 71 = 72 = 73 = yields 



016/2 = V&Ao = y/faPv> = Vfiifo, 
which after scaling results in 

fc = G i a t + a t°t + a X°i + a i a t + a X a t + a t u i 
+ (<?x a t a \ - °l a t a X + vi<7£ a i - a X a t a X) 

+ 013 ( cr X (7 X <7 l + ^t^t^t - a X a X G \ ~ (7 X (T X' J X) 

+ 0n (o-Xo-x °X ~ ^i^t^t ~ a i a X a t + a X a X <j X) 

+ 2444^ . (74) 

The last case, where all four 7, = 0, may be realized 
by setting to zero just three parameters, say 0ig = 0io = 
012 = 0. This enables us to dynamically eliminate one 
more of the bilinear coefficients, say 0g, and to set all 
the third-order coefficients equal to one by scaling. This 
yields a singular canonic form, 

f c = 03" a+o-+ + cr+a+a+ + a+a+a+ + o-+a+a+ 
+ /3 3 a+a+ + 5 CT 3 + a+ + 6 4<4; (75) 

which still depends on three parameters and allows per- 
mutations. All five sZ-tanglemeters obtained from the 
dynamic equations and dependent on three complex pa- 
rameters are depicted in Fig. [21 

IV. ENTANGLEMENT BEYOND QUBITS 

In this section we show how the nilpotent polynomials 
approach may be extended to describe situations more 
general than assemblies of qubits. In particular, we dis- 
cuss in detail the case of qutrits, each qutrit being trans- 
formed by SU(3) or SX(3,C) groups, and construct ap- 
propriate nilpotent polynomials F and / = In F for these 
systems. We define the canonic form F c for F and the 
tanglemeter f c . Next, we generalize this technique to 
qudits d-level elements-qudits. 

Remarkably, the nilpotent polynomials formalism al- 
lows us to make contact with the framework of general- 
ized entanglement, introduced in Refs. [Eg, |5^]. While 
the latter provides a notion of entanglement which relies 
directly on physical observables and, as such, is meaning- 
ful even in the absence of an underlying system partition 
(see also H^), an important special case arises in the sit- 
uations where an element structure is specified, but, due 
to some operational or fundamental constraint, the rank 
r of the algebra of local transformations is smaller than 
d — 1, where d is the element dimension. In this context, 
special attention is devoted to spin-1 systems, namely 
three-level systems restricted to evolve under the action 
of spin operators living in the so(3) = su(2) subalgebra of 
su(3). In particular, we show how to introduce nilpotent 
polynomials for characterizing generalized entanglement 
within a single element of an assembly. In such a case, one 



encounters a new kind of the nilpotent variables whose 
squares do not vanish and only some higher powers do. 
We also consider entanglement among different elements 
of such an assembly. 

We conclude by extending the nilpotent polynomi- 
als formalism to the case of multipartite entanglement 
among groups of elements comprising the assembly, that 
is, to the case where different elements of an assembly 
merge, thereby creating a new assembly with elements of 
higher dimensions. 



A. Qutrits and qudits 

In order to describe entanglement among d-level ele- 
ments of an assembly, one needs to invoke the Lie al- 
gebras of higher rank su(d) and their complex versions 
sl(d, C). The construction of nilpotent polynomials for 
such systems is based on the so-called Cartan-Weyl de- 
composition. We illustrate this for qutrits, d — 3, and 
then generalize to arbitrary d. 



Nilpotent variables for qutrits 



Let us start by reminding some basic facts abo ut group 
theory and Lie algebra representation theory |6|1 l70| . 
The d 2 — 1 generators of the algebra sl(d), the complex- 
ified su(d), may be decomposed into 3 sets: 

(i) a set Ti. of r = d — 1 linearly independent, mutu- 

ally commuting generators of a Cartan subalgebra 
span(H) (r is the rank of the algebra and speoi(TL) 
is the vector space spanned by Ti). In a faithful 
matrix representation, the most natural choice for 
Cartan generators are traceless d x d diagonal ma- 
trices. 

(ii) a set {£} of d(d— 1)/2 "raising" generators spanning 

a nilpotent subalgebra span(£) C su(d). The ele- 
ments of {£} and of sp&n(£) can be represented by 
the matrices with nonzero elements only above the 
main diagonal. 

(iii) d(d — l)/2 Hermitian conjugate "lowering" genera- 
tors spanning a nilpotent subalgebra spanjjF} rep- 
resented by matrices with nonzero elements only 
below the diagonal. 

In the case of su(2), each of the above sets contains 
a single generator: H = {ct 3 = <j z }, £ = {er + } and 
T = For su(3) having eight generators represented 

by the eight Gell-Mann A a matrices [69| , the Cartan sub- 
algebra involves two generators usually chosen as 



I 
A 3 = I -I 




1 



I 

1 I . (76) 
v- 3 \ -2 



23 



The basis £ for the raising nilpotent subalgebra is com- 
prised of 3 elements, 



-t~,u+ = s+. 



In this case, a generic qutrit 



A 1 + jA 2 




A 4 + iX 5 




A 6 + iX 7 




(77) 



(78) 



(79) 



while T includes their Hermitian conjugates. 

The elements of £ and T are the root vectors of su(3). 
It means that for e 6 £, / € T, and h 6 span(H), 
the commutator [e, /i] is proportional to e, and the com- 
mutator [/, h] is proportional to /. The subalgebra 
£ is nilpotent, meaning that the multiple commuta- 
tors [ei, . . . , [e p _i, e p ] . . .] vanish starting at some level 
p — 1. For su(3) with only one nontrivial commutator 
[s + ,t + ] = u + , double commutators already vanish. For 
su(d), they vanish at the level d— 1, coinciding with the 
rank of the algebra r. 

A generic pure state of a qutrit may be represented as 

|V) = Vo|0) + Vi|l) + $j|2) = (Vo + Vi* + + V>2« + ) |0) , 

(80) 

with 



state would be expressed as 

\lf>) = + Vl + ^2M + ) |1) = (ifot- +4>!+ |1> , 

(82) 

The choice of Eq. (|80|) is more natural for discussing 
qutrits, whereas the choice made in Eq. I|82l) allows one 
to associate A 3 with the spin projection operator S z for 
a spin-1 system, with the vacuum state |0) = |1) natu- 
rally corresponding to the lowest eigenvalue —1 of this 
operator. The state |1) is so-called extremal weight state 
for the representation and is discussed in Sect. IVC in 
more details. In any case, this choice is mainly a mat- 
ter of state labeling, which in most cases is dictated by 
convenience and can then be done accordingly for each 
particular physical problem. 

In analogy with Eq. I|13|l . the state of n qutrits can 
then be written as 



l*> = 



E 

,+ 



=<>} 



(83) 



where |0) = |0, . . . , 0) denotes the reference state. Our 
next steps are: (i) to consider the nilpotential / = ln_F; 
(ii) to bring it into the canonic form specified by the 
requirement of maximum population of the state | O) fol- 
lowed by the requirement of maximum population of the 
state |1) = |1 . . . 1), as discussed at the end of Sect. Ill Al 
and (in) to normalize \^c) by the condition (0|^c) = 1- 
The tanglemeter fc thereby obtained provides one with 
simple entanglement characteristics and relevant insights. 
Let us explain how this construction works for an assem- 
bly of two qutrits. 



10) 



2. Entanglement of two quints 



Equation (|80|l generalizes a similar representation for the 
qubit pure states extensively discussed so far. The oper- 
ators t + ,u + belong to £ and commute. Note that even 
for larger d, such a set of d — 1 — r commuting nilpo- 
tent operators always exists [7l| , allowing one to express 
a generic qudit state as a first-order polynomial of com- 
muting nilpotent variables. 

The state |0), which by definition is annihilated by all 
"lowering" generators of J 7 , plays the special role of the 
Fock vacuum or reference state. But the choice Eq. I|81|l is 
not unique. Actually, any qutrit state can serve as a valid 
reference state. For each |6) = f7|0), U G SU(S), the 
analysis above applies if one merely changes (conjugates) 
accordingly the Cartan-Weyl decomposition and defines 



u = uwj- 



£ = U£U- 



T = UTU- 



Again, span(£) and span(F) are nilpotent; again, |6) are 
annihilated by the elements of etc. For example, if 
the choice |0) = |1) is made, the elements of £ are s + — 



We select |0) = |0, 0) as the reference state. The rele- 
vant nilpotent variables are uf and tf , where the index 
i = 1,2 labels the qutrits. An arbitrary assembly state 
can then be written as 



\^)=F(uf,tf)\0) = (Voo + Voi^ 

+ Tptft^uf + ^lU^tf + -022W^Wi")|O) 



(84) 



The requirement of maximum population of |0) dictates 
that the linear terms in F(uf , tf) vanish. This is proved 
by inspecting the variation of the population under local 
unitary transformations, similarly to what was done for 
qubits (see Eq. J!}). Imposing also our normalization 
constraint dictating that the expansion of F starts with 
1, we bring F c into the form 

F c = l+ant^tf + ai 2 t2U^ + a 2 \u^t'l + a 2 2U2ul . (85) 

Thus, we are left with the reference state and four excited 
states. The form in Eq. (|85|l is invariant with respect to 
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the subgroup SU 2 {2) © SU\{2) of local transformations 
which can mix the levels |1) and |2) for each qutrit, but 
preserve the population of the state |0). Using this free- 
dom, we further restrict the canonic form by acting on 
I*) with the transformations of SU 2 {2) ® SUi(2) that 
maximizes the population of the state |11) = t^tf |0), 
while the population of the reference state |0) is al- 
ready maximized. This eliminates the terms oc t^u^ and 
oc u^ti in F c . Thus, we finally obtain 

F c (u+,t+,ut,t+) = 1 + a n t+t+ + a 22 u+u+ . (86) 

The polynomial F c depends only on two complex param- 
eters an, which can be set real and positive by a local 
phase transformation, 

exp (ij 2 X% + i5 2 X% + ijiXl + iSiXf) . (87) 

Thus, two real parameters are sufficient for characteriz- 
ing entanglement between two qutrits, as illustrated in 
Fig. ^ b) . This is consistent with the result obtained 
by a straightforward application of the bipartite Schmidt 
decomposition. 

When counting the number of invariants iVj for a two- 
qutrit state with the help of the expression 

D su = 2 • 3 n - 8n - 2 , (88) 

we find Nj = for n = 2. As for the two-qubit case, this 
number differs from the actual number of independent 
parameters because the phase transformation of Eq. I|87|) 
has more parameters than the number of the coefficients 
in Eq. Q85[) . and some of them act on the coefficients in 
the same way (cf. Eq. (|10|l and the discussion there- 
about). The counting given by Eq. 188(1 holds for n > 3. 

3. Entanglement in a generic audit assembly 

The above analysis suggests the following generaliza- 
tion to an assembly of n qudits. Let di denote the local 
dimension of the i-th element with the associated su{di) 
algebra. For each element i, we choose a reference state 
|0)j and perform the corresponding Cartan-Weyl decom- 
position Tii © £i © J-i for the generators of this algebra, 
in such a way that f\0)i = for all / € Ti). The most 
convenient choice for |0)j is the state with only the low- 
est level occupied. One may choose a basis in the qudit 
Hilbert space involving the vacuum state |0)j and the 
"excited states" , such that each basis state represents a 
joint eigenstate of all generators X^ in the Cartan subal- 
gebra (/tj = 1, . . . , di — 1). The eigenvalues of {A"*} thus 
provide good quantum numbers labeling the qudit state. 

Next, we choose a set of commuting nilpotent genera- 
tors {^i'} C £j that may be employed as nilpotent vari- 
ables. To be specific, let us enumerate these variables 



such that = v\ corresponds to the first excited 
state of i-th. element, |2) i = vf |0) i to the second, etc, A 
polynomial F {V" 4 } of these nilpotent arguments charac- 
terizes a generic quantum state of the assembly as the 
latter may be obtained acting by the operator F on the 
corresponding assembly reference state |0) = Yii 

In analogy with the procedure for two qutrits described 
above, we can, using local SU(di) operations, maximize 
the population of the vacuum state and eliminate the 
terms linear in the nilpotent operators. The function F 
acquires the form generalizing Eq. I)85[l. 

= ! + v i v i +^ijk v i v i v h ---y ( 89 ) 

where repeated indices imply summation. The state is 
normalized to unit amplitude of the reference state, as 
earlier. As a next step, we maximize the population of 
the symmetric state |1, . . . , 1) using the transformations 
of the subgroup ®iSU{di — 1), where the tensor prod- 
uct is taken over all i with di > 2, which preserve the 
reference state. At the third step, we maximize the pop- 
ulation of the state |{fci}), where ki = 2 for the elements 
with di > 2 while for qubits the label ki is "frozen" at 
the value ki = 1, and maximization is done using the 
transformations of the subgroup ®iSU(di — 2) that af- 
fect neither the reference state nor the first excited state, 
and the tensor product involves now the elements with 
di > 3. If the assembly involves 5-level elements, we 
are allowed at the next step to maximize the population 
of the state where fc, = 1 for qubits, ki — 2 for 

qutrits and ki — 3 for the elements with di > 3, etc. 

For example, for an assembly of a five-level system, a 
four-level system, and a qutrit, one should consecutively 
maximize: 

(i) the population of the state |0, 0, 0) by the transfor- 
mations from SU(5) ® SU(4) <g> SU(3); 

(ii) the population of the state |1, 1, 1) by the trans- 
formations from SU(4) ® SU{3) <g> SU(2); 

(Hi) the population of the state |2, 2,2) by the trans- 
formations from SU(3) ® SU(2); 

(iv) the population of the state |3, 3,2) by the trans- 
formations from the remaining SU(2) mixing the 3-d and 
the 4-th excited states of the five-level system. 

This is all illustrated in Fig.0] 

This procedure eventually reduces the nilpotent poly- 
nomial F to the canonic form of Eq. I(89|l . where some 
coefficients a^ Kj , a K ^ ,K ' k now vanish. In order to better 

l J IJK 

understand the pattern, consider two examples: an as- 
sembly of three qutrits, and of two qutrits and a qubit. 
In the first case, the canonic form is 
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F c = 1 + aont^ti + aout^ul + aoaiw^i" + 0022^2 "i" + a ioi*3~*i" + 0:102*3" wj" + 0201^3^ + 0202^3 + ono*3~*2~ 

+ 01120*3^2 + a 2WU^t2 + a 2 20ut U2 + aillfrj" t^ti + Ce^t^U^U^ + a 2 l2U^t2U^ + a22lU^U2tf + a 2 2 2 u^ u^u^ l 

(90) 

whereas for two qutrits and a qubit (labeled by the index 3), it looks as follows 

F c = 1 + otit^ti + 0:5*2 u i + a 7 u 2 H + ct^u^u^ + aio*^"*^" + ant^iii + 012*^"*^ + Ois^uJ + 013*3 t^n + 0:17*3 U^u^ , 

~ " [911 



r 



with the identification v\ = tf and v\ = uf . In Eq. (|90|1 . 
the indices of a explicitly label the individual states of 
the qutrits, while in Eq. I|91[l the notation relies on the 
decimal representation of the base-3 numbers associated 
with these indices, 112 — > 14, etc. The canonic forms are 
defined, as before, up to phase factors of the nilpotent 
variables . 

The forms of Eqs. (|90I91|) do not involve linear 
terms. Neither do they involve any term proportional 
and (for three qutrits) u^t^tf. In 



to tft^uf, 



4° 



other words, the amplitudes of the states |1, 1, 2), |1, 2, 1) 
and, for three qutrits, |2, 1,1) vanish. This vanish- 
ing is achieved at the second stage of our procedure 
maximizing the population of the state |1,1,1) by the 
transformations from ®i =1 SU(2)i for three qutrits, and 
from SU 2 (2) <S> SUi(2) for two qutrits and a qubit. In- 
deed, from the viewpoint of the remaining SU (2) trans- 
formations, the state |1) may be regarded as a ref- 
erence state, and amplitudes like 1 1,1,2) vanish for 
the same reason why the amplitudes |1,0, ...,0) van- 
ished when the population of \0) were maximized. It 
is clear that, e.g., for an assembly of three four-level 
systems, where the procedure involves 3 steps, the 
canonic state has vanishing coefficients of the basis states 
|2,1,1),|1,2,1),|1,1,2),|3,1,1),|1,3,1),|1,1,3),|3,2,2), 
|2, 3, 2) and |2, 2, 3). Generally, if the procedure involves 
the maximization of the population of the state |{fcj}), 
the vanishing amplitudes are |{^}), where the sets {ki} 
and {k'j} differ only in one position i$ and k[ Q > fcj . 

Going back to Eq. i|9"U|) . we observe that the canonic 
polynomial depends in this case on 17 complex param- 
eters, 6 of which can be set real by a proper choice of 
the phase factors of the nilpotent variables uf and tf. 
This coincides with the number Nj = 28 suggested by 
the asymptotic formula given in Eq. I|88|) valid for n > 2. 

The coefficients a^"^ 3 can be treated as invariants 
characterizing the qudit entanglement. But the coeffi- 
cients of the tanglemeter f c = hiF c provide, as we al- 
ready saw, a more direct and physical description. In 
particular, the criterion of Eq. (|18|) . indicating whether 
two groups, A and B, of a bipartition are entangled, may 
be straightforwardly generalized. We have 



= 0, i e A, m e B . 



As before, this criterion holds even for a noncanonic 



nilpotential. 

The construction of the sZ-tanglemeter may be accom- 
plished by analogy to the qubit case. One can derive 
and directly employ the dynamic equations for /, elimi- 
nate the states adjacent in the sense explained above to 
\k, . . . ,k) by a proper choice of the transformation pa- 
rameters P/% thereby obtaining the sZ-canonic form: 



(92) 



Like it was the case for the su-tanglemeter, many of the 



coefficients ft- 



Actu- 



Pijk J 1 etc vanish, however, 
ally, as SL transformations have more parameters than 
iSCZ-transformations, we can now bring to zero more co- 
efficients than in the SU case. Namely, the sum on the 
right hand side of Eq. (|92|l contains no terms correspond- 
ing to the states directly coupled to \k,...,k) by any 
single local transformation, and not only those involving 
"higher" states with k! ^ k (for simplicity, we discuss 
here only the case when all elements have the same di- 
mension). Without explicitly presenting the correspond- 
ing dynamic equations, we illustrate this idea in Fig. 0] 
for the case of qutrits. As one can see there, the generic 
SL(3, C) canonic form of / contains 8 complex parame- 
ters, 6 of which may be specified by a proper choice of 
the scaling factors. Thus, 

f C = Pg {ujt+ + U+t+ + t+U+) + t+t+tt 
+ 0u (*2~ u l~ + u 3~*l~ + *3~ u 2~) + u t u 2 U t ■ 



B. Generalized entanglement and generating 
functions 



We now consider a situation where the local operations 
available for qudits are restricted by some operational or 
fundamental constraint, so that they cannot ensure uni- 
versal local transformations. In particular, we focus on 
the case where the restricted local operations form a sub- 
group ST/ (to) C SU(d), with to < d, of the full unitary 
transformation set. As mentioned, this is a special rele- 
vant instance of a more general, subsystem-independent 
entanglement setting (generalized entanglement) formal- 
ized in Refs. [5^, |59| ■ One of the main implications of 
the latter approach is, in turn, to point to an intimate 
connection between entanglement and so-called general- 



20 




FIG. 4: a)A possible strategy for identifying the tangleme- 
ter f c for three qutrits. The choice of the parameters P* of 
local unitary transformations is done in such a way that the 
transformations sequentially eliminate the terms in F c that 
are coupled to the states \k, fc, fc) by just one local operator. 
The SU(3) transformations maximize the population of the 
vacuum |0) = 1 0, 0, 0) . Next, the SU(2) transformations not 
acting on the states |0j) shown by dotted circles, and only 
the states shown by dashed and solid circle are mixed, max- 
imize the population of the state |1, 1, 1). The population of 
the state |2,2, 2} cannot be further maximized by local uni- 
tary operations. However, SL(3) and SL(2) dynamic trans- 
formations can eliminate some of the coefficients /3 directly 
connected to such states, b) The same diagrams for the case 
of different elements. 



ized coherent states (GCSs) |72l |73| . which is also useful 
to place our current analysis in a broader context. 

Generalized coherent states may be constructed for 
quantum systems described by a dynamic-symmetry Lie 
group, which typically is assumed to be reductive or 
semisimple: a family of GCSs may be thought of as 
an orbit resulting from application of the Lie group to 
a reference state, which is identified with an extremal- 
weight state in the irreducible representation of the 
underlying dynamical Lie algebra. Within this gen- 
eral setting, the canonical coherent states of a quan- 
tum harmonic-oscillator [74| may be recovered as result- 



ing from the application of the displacement operator 
[exp (az* — a<z)] to the ground state |0). Here, as usual, 
a and denote annihilation and creation operators, re- 
spectively, and z = x + ip G C stands for displacement 
in the phase space (x,p). The relevant dynamic groups 
are generated by the Heisenberg-Weyl algebras, cither 
f)3 = {I, a, a'} that allows for arbitrary displacements, 
or h 4 = {I, a, a\ a' a} that in addition involves the phase 
space rotation exp (iipa^a) by an angle tp |75j |. Physi- 
cally, the most representative manifestation of the co- 
herent states is the fact that external interventions re- 
alized via a] and a cannot change the structure of the 
wave function expressed in terms of the Wigner function, 
but just displace and rotate it in the phase space. One 
can also include the operation [expi (s(a^) 2 + s*a 2 )] of 
squeezing seCby remaining within a finite-dimensional 
(so-called "two-photon") dynamic algebra. Clearly, the 
algebra describing arbitrary transformation on the oscil- 
lator Hilbert space is infinite-dimensional, and destroys 
the manifestations typical of coherent states. 

Interestingly, the mathematical techniques employed 
to describe GCSs are extremely close to the ones we use 
for nilpotent polynomials: any GCS is constructed from 
a minimum-weight reference state of a semisimple Lie 
algebra by application of a function of the appropriate 
Cartan-Weyl raising operators which, for a finite dimen- 
sional Hilbert space, are indeed nilpotent. One can al- 
ways construct GCSs of a qudit, relative to its full su(d) 
algebra of transformations. The resulting family is, of 
course, identical with the set of possible pure states of 
the system. For a qudit assembly with n elements, it 
was shown in |58l | that conventional wnentangled states 
are identical with the set of GCSs under arbitrary lo- 
cal transformations, that is, with respect to the algebra 
®" =1 su(e?j), with rank R — J2i(di ~ !)■ Geometrically, 
such GCSs are extremal (in the convex sense) in the set 
of states which may specified only through expectation 
values of arbitrary local observables. However, such ex- 
tremality property may no longer be fulfilled upon re- 
stricting the transformations on each element to a proper 
subalgebra, mi < di - implying the possibility of general- 
ized entanglement relative to the restricted observable set. 
Here, we further develop the connection between entan- 
glement and GCSs in terms of the nilpotent polynomials 
approach. In a way, we can say that we consider en- 
tanglement related to GCSs both within each element as 
well as entanglement among GCSs of different elements 
of the assembly, and construct in each case appropriate 
characteristics based on nilpotent variables. 

The first step toward accomplishing the above goal is 
to obtain a proper description of the group of restricted 
local transformations, by embedding it as a subgroup into 
the full-rank = di — 1 group of local transformations. In 
other words, among the generators of full group we have 
to specify the linear combinations that correspond to the 
generators of the restricted local transformations. To this 
end, it turns out that so-called generating function tech- 
nique offers a convenient tool. Let \^>) be a state of the 
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assembly normalized to the unit amplitude of the vacuum 
state, and apply the operator exp J2i X T { V T^ ; where 
the nilpotent operators v*' correspond to the Cartan- 
Weyl decomposition of the full su{di) algebra of i-th el- 
ement, and x** € C. We finally project the result onto 
the vacuum state and obtain 



F(K'})HO|exp|V*H^) T I*) 



\t 



(93) 



which is the generating function for the set of variables 
{x^}. This coincides with the function F ({v^}) given 
by the nilpotent polynomial F where all the operators 
o~ + are replaced by the corresponding variables x. 

As long as the set U** corresponds to the full algebra, 
introducing the generating function adds nothing new 
to the characterization of entanglement. The situation 
changes when we consider subalgebras of rank < r, . 
The generating function takes then the form 



F(K'}) = (0|exp[] 



I*) 



(94) 



where the commuting nilpotent operators {/J,** } belong to 
an algebra of rank m^. The elements ^ are linear com- 
binations {^i*} C span-fi/^} of the nilpotent elements in 
the Cartan-Weyl decomposition L z L + L~. Among 
the Cartan generators {A K }, we also have to single out 
a subset {v K } C L z of the operators corresponding to 
the choice of {^j- Note that linear combinations [if* 
of the nilpotent operators v?* are also nilpotent, how- 
ever of higher order, that is, such that (^^' , ) Pi = for 
some integer p** > 1. For example, for /i = t + + s + of 



Eqs. CZJ-CSl, one finds [i 2 



^ 0, /r* = 0. Therefore, 



unlike all cases discussed so far, the generating function 
F {{ X T}) ma V contain some higher powers of the vari- 
ables x^' , along with the terms linear in x^' . 



The generating function F ({x^ }) may also be reduced 
to canonic form. This has to be accomplished, however, 
only by resorting to local unitary operations which belong 
to the subalgebra of restricted local transformations, 



F c ({<*}) = <0| exp 



exp [i5»r + r? + {rtt w 



(95) 



r 



Note that one may also generalize this approach to the 
case of restricted SL transformations, where in place of 
Z"*, R1\ and (-Rf*)*, one has to substitute in Eq. l|9^|) 
independent complex parameters. The corresponding 
nilpotent polynomial F c ({fJ-^ 1 }) may now contain powers 
of the nilpotent variables, which is a signature of a very 
important fact: Entanglement is no longer necessarily as- 
sociated with different subsystems, but may occur within 
each single element as well. 



C. Examples: Generalized entanglement for one 
and two spin-1 systems 

We first consider generalized entanglement in the sim- 
plest example of a three-level system (see Refs. [l9l 158^ 
and show how this result may be interpreted in terms of 
nilpotent variables. Though a three-level system corre- 
sponds to a full su(3) algebra, we consider it here as a 
spin-1 system that is, concentrate on a situation where 
the physical observables are restricted to the subalgebra 
su{2) of spin operators, S± = S x ±iS y and S z . Note that 
the latter are equivalent to the u ± +t T , and A 3 generators 
of su(3), respectively. The spin states are characterized 
by the eigenvalues of v = A 3 , which is the only Car- 
tan generator of the su(2) subalgebra. The spin-down, 
lowest-weight state | — 1) is chosen as the reference state. 
The operators s + and t~ form the commuting nilpotent 
subalgebra {v} of su(3) and give two other "excited" 



states |0) = t~ | — 1) and |1) = s + | — 1). They character- 
ize the quantum states of the three-level system according 
to the relation |*) = (1 + a s s + + a t t~) |-1). The oper- 
ator S+ — u + + t~ is the only element of the nilpotent 
{/x"*}-subalgebra su(2) c su(3). 

Now we show that the state |0) is generalized entan- 
gled with respect to SU(2). Indeed, by the unitary ma- 
trix ^-™s^l\f2 gtatg |q) may be transformed to the 
state with the maximum vacuum population \ipc) = 
(| — 1) + |1)) /V2, which is evidently different from the 
reference state | — 1). The corresponding canonic state 
normalized to unit reference state amplitude reads 



|* c ) = (l + s+)|-l) =exps+|-l) 



(96) 



hence f c = s + and F c = 1 + s + . We now construct the 
generating function of Eq. (|93|l by employing S+ = u + + 
t~ as the only element of the nilpotent }-subalgebra 
of su{2) for our system. This yields 



F(x) = (-l|exp x (u + +t~) 

= (-1| exp [x (u~ + t + ) 
= l+x 2 /2 . 



exp s | — 1) 



exp s 



l-D 



(97) 



The presence of the quadratic term x 2 /2 is the signature 
of generalized entanglement. 

One way to understand the meaning of this "self- 
entanglement" in the state |0) is to see it as a consequence 
of the fact that the operators S±. z of su(2) cannot lift 
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the degeneracy of the two eigenstates of the operator As, 
which together with A3 = S z labels the states in the unre- 
stricted su(3) algebra. In other words, within the group 
of restricted local transformations, the transition from 
the state |0) can only access the state (| — 1) + |1)) /V2 
but not (| — 1) — |1)) / \/2, hence the amplitudes of the ref- 
erence state | — 1) and that of the state |1) are fully cor- 
related. Alternatively, one can say that no SU(2) trans- 
formation is able to connect the state | — 1), which is a 
SU(2)-GCS and is unentangled, to the state |0) which is 
a a SU(3)-GCS but not a 5[/(2)-GCS. 

For a generic SU(3) state (1 + a s s + + a t t~) \— 1), the 
nilpotent SU(2) polynomial 



F 



Fix) 



is given by the generating function 
F(x) = {— 1| exp x(u + +t~) exp(a s s" 

„2 



1 + atX + 



a s x 



(98) 

a t t~)\-l) 
(99) 



The corresponding canonic nilpotent polynomial of 
Eq. I|95|) and the tanglemeter take the form 



F C (S+) = 1 + a' s S 2 + /2 
f c (S + ) = a' s S\/2 , 



(100) 



respectively, where the phase of a' s can be set to 0. There- 
fore, generalized entanglement of a single spin-1 is char- 
acterized by a single real parameter. 

Consider now a second example, that is, entanglement 
between two spin-1 |19L |58j, which we describe as two 
three-level systems subject to the action of SU(2)(BSU (2) 
local operations. By analogy to the single spin-1 case, we 
chose nilpotent variables = {5i )2 } in stt(2) © su{2) , 

where 5i l2 = (5+) 1,2 = «^ 2 +*i> and = {O^)^}- 
The state | — 1, — 1) is chosen as the SU (2) reference state. 
As before, the operators s x , t x , sf 2 , and t X2 are the 
nilpotent variables in the full su{3) © su(3) algebra. 
A generic quantum state for a two-qutrit assembly 



|^) = (1 + a S ;i s t + a t;it\ + a s - 2 s 2 + a t -2t 2 + a s , s s 2 s^ + a t ^s 2 ^t 1 + a t , s t 2 + a t>t t 2 t 1 ) |-1, -1) , 
is now characterized by the nilpotent polynomials 

F ({s + , t~}) = 1 + a s . A sf + attjitf + a s - 2 s 2 + a t -2t 2 + a StS s^sf + a t}S s^t x + a tjS ^s+ + a t , t t 2 t x , 



f({s + ,t }) = a s;1 sf + a t ;xtx + a s -2 s 2 + a tah 



+ (a s , s - a s; 2a s; i) s 2 s+ + (a t , s - a s . 2 a t .i) s%t x 



+ {a s ,t - a t -20L s -i) t 2 s x + (a t ,t - a t . 2 a t -i) t 2 t 1 ) , 
whereas for the SU{2) characterization, Eq. I|94|) gives the generating function 



F(x,y) = (-l,-l|exp x(u+ + t x ) 



t -\t 



y{ut + t 2 -) f \F({s+,t-}) 1-1,-1) 
I 



(101) 

(102) 
(103) 

(104) 



Direct calculation yields 

2 2 
x y 

F(x,y) = l+a s; iy + a t ,ix + a s , 2 — + a f;2 y + a t . t xy 

2 2 2 2 

+ a t , s x— + a st —y + a s , s — — ■ (105) 
This finally results in the nilpotent polynomial 



F(S 1 ,S 2 ) = l + a s . 1 - 



Ott;lSl 



a S ;2- 



si 



and | — 1,0) vanish, thus one obtains 

o 2 q2 q2 c2 

F c (S lt 5 a ) = 1 + a s:1 ^- + a s , 2 ^ + ^,^2 ( 107 ) 



02 

1 2 q 



otstS: 



5? 



f c (S 1 ,S 2 ) = P s ,i^ + 0t,tSiS 2 + (3 st ^S 2 



02 



Si 



■A 



' 2 2 



C2 c 2 c2 o2 

, °2 D l . °2 a 1 c 1 1 a a 
+ Os.syy + a t,s~2 01 + a st b 2 — + a t ,to 2 oi , 

(106) 

where the subscripts + of the nilpotent variables S are 
implicit. In the canonic form maximizing population of 
the reference state, the population of the states 1 0, — 1) 



where f3 s>s = a SyS ~2at,t 2 , and all other (3 — a. As before, 
by exploiting the freedom of phase transformations on the 
nilpotent variables, the parameters a s -\ and a s;2 (or f3 s -i 
and /3 S;2 ) characterizing generalized entanglement within 
each of the three-level systems can be set real, and we are 
left with four complex numbers a s . s , at, s , a s t, and at,t 
characterizing the generalized inter-spin entanglement. 



29 



One may also ask the following question: How can we 
characterize generalized entanglement under SL trans- 
formations? For qubits, the resulting classification is 
based on the dynamic equation (|54|l and the conditions 
Eq. (|65[) - l|66l) of the exponential decrease of the chosen 
coefficients in the course of controlled local SL(2, C) dy- 
namics. In order to suggest a strategy for characterizing 
an ensemble of spin-1 elements, we note that an analog 
of Eq. H54fl may be derived for nilpotent polynomials on 
Si with the help of the analogs of Eqs. (|41|) . A require- 
ment of eliminating certain terms in / by a proper choice 
of the local SL(2,<C) control parameters can also be im- 
posed by analogy. In particular, the part of /, which 
contains only multi-linear terms in Si may be reduced 
to the SL(2) canonic form which, up to the replacement 
&i~ i ► Si 7 is identical to that of qubits. This specifies the 
canonic form of the entire /, whereas the terms contain- 
ing at least one Sf factor comprise a generic polynomial 
with the coefficient specified in the process of reducing 
the multi-linear part to the canonic form. 

D. Transformation of the nilpotential under change 
of partition 

We next consider a situation which, in a sense, is the 
opposite to the above-discussed scenario of generalized 
entanglement, whereby the ranks of the algebras em- 
ployed for the entanglement classification exceed the di- 
mensions of the local Hilbert space of the elements. This 
is the case of a partition of a composite system, where 



each part is composite by itself and may contains mul- 
tiple elements. In other words, from an initial assembly 
we compose a new one by considering groups of elements 
as new elements and by describing the quantum state of 
each group by a single, collective quantum number. One 
can say that the new assembly results from merging of 
elements in the old one. 

We begin with the example of n-qubit system charac- 
terized by the nilpotential of Eq. I|15|) , 

n 

f(Wt})= £ a^.^IIK)"' ( 108 ) 

{fci}=0,l i=l 

partitioned in three parts A, B, and C, each of which con- 
tains tla, Ti b i and nc qubits, respectively. The particular 
case C — recovers the bipartite setting. The new as- 
sembly thus consists of three elements with d,A = 2 nA , 
d B = 2 nB , and d c = 2 nc . 

By exploiting a standard Hubbard-Stratonovich pro- 
cedure |76j| . we can represent the polynomial F = e* 
corresponding to Eq. (|108|l in the form of an integral 

F({af})= /e^»fte-l-l 2+2 -^ d " Z ' d " < , (109) 
J ti 

where the integration has to be performed independently 
over both the complex variables Zi and their complex 
conjugates z* . This suggests a straightforward separation 
of the system into three parts, 



F |0) = / exp 



J2 ztaf |0) c ex P J2 <°t l°>B ex P £*rt I0^e/««'>>ne-"<' a ^^, 

■ A+ng + l ) [i=n A +l J U=0 J i=l 



where \0) a-b-c denote vacuum states of the new elements, 
which we still can choose as product states of qubits in- 
cluded within the new elements. 

Let the Cartan subalgebras and the commuting nilpo- 
tent elements of the su(2 nA ), su(2 nB ), and su(2 nc ) al- 



gebras be denoted by {\ K /} {A^ B } {v B B }, and 

{A£f } {v£ C }, respectively, and let the state of the as- 
sembly be characterized by a nilpotent polynomial on v K 
upon noticing that (V K ) 2 = 0. This yields 



l + E<°lc(^ C ) f e 



En ^ 
... . i=n A +n B + l z i 



\Q)c"c 



kC 



IT 

i=l 



"Nil 



2 d™Z,-d n Z* 



(111) 
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One may rewrite this in the form 



2"A ; 2"B ;2"C 

F, r KA «B KC 1 \ \ ^ 



(112) 



where the expression 



/" f " A . 1 f TI A+U B ~\ 

a KAlKB> na = / (0U W) f ex P E l°>A <°U ex P E *?"< + 

U=0 J U=7iA+l J 



IC 



(^) f exp £ <<r+l|0) c e^««*»ne- 



ni- 

2 d"z 2 d"z* 



i=l 



7T 



(113) 



r 



explicitly gives the coefficients of the new nilpotent poly- 
nomials characterizing entanglement in the new assem- 
bly. Here, the identity operators v\, v^, and Vq are in- 
cluded in the sets , {^s B }, and {z^°}, respectively. 
Note that for bipartite and tripartite entanglement, the 
expressions for a and for f3 are identical, provided that 
one eliminates terms of F linear in v by local transforma- 
tions SU{2 nA )®SU{2 nB )®SU{2 nc ), and normalizes the 
reference state population to unity. Also note that for 
a bipartite case, the sum p KAtK > A = S« B 2 KA , KB 5*^ refl 
gives the density matrix of part A normalized to unit 
population of the reference state. 

The generalization of expressions Eq. I|112|) - (|113|l to a 
larger number of new elements is straightforward: 

2™A;...;2"W 
K A , . . . ,K \y —0 

2 nA ;...;2 nw 

/(KV ■ ■ » Kn) = E ~^A,.,^y R A A 



w 
W > 



where the integrand of Eq. (|113fl for a K 
tains more factors 



(114) 



now con- 



ic 



^) t exp{E i ^ 



■}|o> 



A" 



with i running over all qubits included in the new element 
set K. Once the reference state amplitude is normalized 
to 1, the nilpotential / may be found by direct calculation 
of lnF. This yields relations among a and /3 that even 
for the canonic state do not coincide when the assembly 
comprises more than three new elements. 

One can finally raise the following question: What 
could be a reasonable choice of a state marking an or- 
bit in the case where the dimensions of elements be- 
come large, and the exponential complexity makes the 
identification of a global maximum of the reference state 
population an intractable problem? One of the possi- 
bilities is to rely on dynamics and manipulations exclu- 
sively with the nilpotential /, aiming at the elimination 
of the coefficients /3 of all the states connected to the 



states \k, . . . , k) by a single local transformation, as illus- 
trated in Fig. 21 However, we note that such a choice, 
though providing a unique characterization of entangle- 
ment, may yield canonic forms not corresponding to the 
maximum population of the reference state. Hence it 
can be ambiguous, leading to different polynomials for 
the same orbit, as it is the case of the example presented 
in Eq. (25) of Ref. j^. Still, such a choice might have 
some advantages in view of "operational compatibility" 
important for large ensembles, since it relies exclusively 
on the manipulations with the extensive polynomial / 
and does not invoke F, thereby avoiding the need for the 
exponentially long procedure of F <^> f conversion. An- 
other advantage is that the generalization to the SL case 
is straightforward. 



V. SUMMARY AND OUTLOOK 

We conclude by summarizing the general idea of our 
entanglement description in composite quantum systems, 
based on the selection of a product reference state and 
on the introduction of appropriate local nilpotent opera- 
tors. We also discuss the results obtained with the help 
of the nilpotent polynomials and mention possible future 
applications of this technique to other physical problems 
relevant to quantum entanglement and quantum infor- 
mation science. 

Dealing with a physical system composed of distin- 
guishable parts, we need to specify among which subsets 
of the parts we wish to consider. In order to avoid confu- 
sion, we consider each of these subsets as a single element 
and characterize it by a single quantum number. We call 
an assemply, the collection of all the elements. 

The main purpose for introducing the nilpotent vari- 
ables technique is to obtain a characterization of entan- 
glement by extensive quantities that are sums of the 
characteristics of the unentangled parts of the system. 
This approach critically relies on an important property 
of nilpotent polynomials: any analytical function of a 
polynomial is also a polynomial, whereas a key role is 
played by the logarithm function, which enables one to 
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relate the product state to a sum of independent terms, 
each of which represents a part of the system unentanglcd 
with the rest. The extensive characteristic, nilpotential 
/, that emerges, for a product quantum state is indeed 
a sum of nilpotentials of the unentangled parts, whereas 
the presence of entanglement among the elements of dif- 
ferent parts is represented by the corresponding cross- 
terms. Verification whether or not such terms are present 
in the nilpotential thus serves as the entanglement crite- 
rion. 

By the very meaning of entanglement among elements 
of an assembly, these characteristics should be insensi- 
tive to local transformations, which result from arbitrary 
reversible (both unitary and not unitary) individual ma- 
nipulations of the elements. Therefore, the entire orbit 
of states, that is the manifold of all states of the as- 
sembly that can be reached from a given initial state by 
local operations, should correspond to the same param- 
eter values characterizing entanglement. The number of 
these parameters, that is the dimension D of the orbit 
coset, depends on the type of local transformations al- 
lowed. For both unitary SU and non-unitary SL trans- 
formations, we propose a canonic form of the assembly 
state, which depends on D parameters and serves as the 
orbit marker. For qubits, the canonic form relies on the 
choice of a reference product state of the assembly, the 
vacuum, and on maximization of the population of this 
state via local transformations. 

Given the canonic state of an assembly of qubits, we 
represent it in the form of a polynomial Fc,c({cf}) 
on the raising operators af acting on the vacuum 
state. After normalizing to unit vacuum state ampli- 
tude, we calculate another polynomial, the tanglemeter 
fc,c({o't}) — m ^c : c({o'j + }), which depends on exactly 
D parameters and contains all the information about en- 
tanglement in the form of the coefficients standing in 
front of different nilpotcnt monoms J[ i erf . Subscripts 
C and c show that we remain within SU or SL group of 
local transformations, respectively. In contrast to state 
given by the polynomial F c c, the tanglemeter / Cj c is an 
extensive quantity. In statistical physics, extensive quan- 
tities scale linearly with the size of the ensemble, like 
free energy given by the logarithm of the partition func- 
tion. The tanglemeter has similar additive properties: 
for a composite system representing a set of unentangled 
parts, the tanglemeter equals the sum of the tangleme- 
ters of the constituent parts. Straightforward inspection 
of the second derivatives d 2 f Cy c / dcr^ eA da^ eB allows one 
to check whether or not groups A and B of elements are 
entangled. 

We have presented several examples of tanglemeters 
for systems of a few qubits. Still there may exist a cer- 
tain number of states, for which the tanglemeters, though 
dependent on the same number D of parameters, cannot 
have the chosen structure. These states comprise singu- 
lar classes. In the four-qubit example, we have shown 
in detail how a classification can be constructed by con- 
sidering local infinitesimal transformations sequentially 



diminishing the coefficients of the monoms to be elimi- 
nated. This procedure is described by a set of differential 
equations for the coefficients of / and can be seen as con- 
trol process, when we impose feedback conditions on the 
parameters of continuous local transformations. By prop- 
erly adjusting these parameters to current values of the 
coefficients of the nilpotential, we rapidly arrive at the 
tanglemeter f c in the limit t — > oo. The latter thus ap- 
pears as a stationary stable point of the set of equations 
with feedback. However, the proper choice of the param- 
eters required for elimination of the monoms cannot be 
made in certain domains in the space of the nilpotential 
coefficients, where this procedure fails as a result of the 
"loss of complete controllability" . In these domains, the 
determinants of matrices relating the time derivatives of 
the coefficient with the parameters of the local transfor- 
mations vanish. The cases where one, two, three, or all 
four eigenvalues of the determinant equal zero correspond 
to different entanglement classes. More systematic explo- 
ration of this classification in larger assemblies of qubits 
looks like an interesting prospective task for immediate 
future research. 

We have demonstrated that the nilpotent polynomial 
technique, initially developed for qubits, can be extended 
to d-level elements, qudits, each of which has different di- 
mcntion di. In this case, a larger number of the nilpotent 
variables per subsystem are required for the construction 
of the polynomials. Let us spell out the salient features of 
the proposed description of entanglement for the general 
case. 

(i) We start by choosing a reference state |0)j for an in- 
dividual clement and a corresponding state for the 
composite quantum system, an assemply consist- 
ing of d-level systems, | O) = |0, . . . , 0). All other 
states are obtained by the action on | O) by polyno- 
mials of nilpotcnt local operators. For qubits, these 
nilpotcnt operators are simply the spin raising op- 
erators o~l . For qudits, we invoke d — 1 nilpotent 
commuting operators i/"* that create rf, — 1 excited 
states = v^ 1 |0)i out of the vacuum state |0)j. 
This choice of basis is natural in the framework of 
the Cartan-Weyl decomposition of the su(d) alge- 
bra. Thereby any state of the composite system 
may be represented by a polynomial F{y^ i ) acting 
on the reference state. 

(ii) By applying local transformations to a given assem- 
bly state we bring it into a certain canonic form. 
For systems of qubits, the latter is characterized 
by maximum population of the reference state | O) . 
In the case of qudits, we need to additionally max- 
imize the population of the maximum symmetric 
excited states. When the dimensions di of all ele- 
ments are equal, these are the states |1, . . . , 1), 
\k,...,k)..., \d — 1, . . . , d — 1), while for different 
di, in this sequence,we replace k by di — 1 for the 
elements of the dimcntion di < k. 
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(iii) Any analytical function of a polynomial of nilpo- 
tent arguments is also a polynomial. The log- 
arithm of the canonic polynomial F c , tanglemeter, 
has coefficients that by construction are invariant 
with respect to local transformations. This offers 
a simple way to systematically list all the invariant 
characteristics of entanglement for a system of an 
arbitrary number of qudits. 

Another extension of the technique applies to a case 
relevant to generalized entanglement, where the algebra 
of local transformations has a rank strictly less than that 
suggested by the dimension of the elements. In such 
a situation, the number of nilpotent variables per ele- 
ment is less than d — 1, and moreover, the form of the 
nilpotent polynomial is sensitive to the choice of refer- 
ence state, thus being "reference" and "observable" de- 
pendent. Identification of classes for large elements and 
for generalized entanglement based on the dynamic equa- 
tions are two immediate open questions concerning the 
proposed technique. It would be interesting to identify 
the classes in the case of generalized entanglement and 
to derive the dynamic equations describing the evolution 
of the corresponding nilpotential. This constitutes yet 
another direction for future studies. 

Manipulation of quantum assemblies and, in particu- 
lar, quantum computation, implies application of both lo- 
cal and nonlocal two-particle gate transformations. This 
is also the case for system composed of naturally inter- 
acting elements, such as spin chains, cold Rydberg gases, 
and arrays of cold two-level atoms trapped in a standing 
electromagnetic wave. The Schrodinger equation yields a 
dynamic equation for the nilpotent polynomials - a lin- 
ear one for the polynomial F, which is in turn linear in 
the state amplitudes, and a nonlinear one for / = In F. 
For an ensemble of n qubits, both dynamic equations 
involve ~ 2™ variables, the polynomial coefficients and 
their time derivatives. At first glance, linear equations 
are always simpler. However, this is not necessarily the 
case: we have shown an example of universal evolution 
that yields an equation for /, which is equivalent to the 
classical Hamilton- Jacobi equation in just n dimensions. 
How far does this analogy go? What are the consequence 
for entanglement dynamics and, in particular, for quan- 
tum algorithms? These are intriguing questions to be 
further addressed. 

Finally, one could also conceive applications of the 
nilpotent polynomial technique for analytical investiga- 
tions of entanglement in correlated many-body systems - 
including the problem of better characterizing quantum 
phase transitions [l9| - and extend the present consider- 
ation to the case where non-unitary decoherence and dis- 
sipation effect are included in the dynamics. In an even 
broader context, one can think of employing this tech- 
nique as a tool for establishing relations between entan- 
glement dynamics ^^,0] and the known exactly solvable 
problems of statistical mechanics and complex quantum 
systems [77tl78| . 
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VI. APPENDICES 

A. Dimension of cosets 

a. Two qubits The fact that the expression Eq. (0J 
is invariant implies that the counting 2 n+1 — 3n — 2 of the 
number of invariants breaks for n = 2. That means that 
out of 3n = 6 local transformations, only 5 act faith- 
fully (non-trivially) . In other words, there is a certain 
local transformation which leaves a generic wave func- 
tion intact. On can amuse oneself and find it explicitly. 
A generic infinitesimal local transformation is 

Leaving only the components V'oo and ( as dictated by 
the canonic form), and writing P^c^ = P* a z + P~a + + 
P^p~ with r = 1, 2, we arrive at the system of equations 

-i<fy»ii = {PI + P^n = 
-i6iP 00 = -(P 1 z +P2)ipoo=0 

-iSipoi = Pf Voo + P^^n = 

-iSfao = P^oo + PMi = (116) 

From the last two equations, it follows that Pf = P 2 + = 
(for generic ipoo , ^n), while the first and the second ones 
give the same condition Pf + P| = (or fa + fa + <p i + 
if2 = in the notation of Eq. I|1U|)1. i.e. we have, indeed, 
a one-parametric set of local transformations that leave 
ipij intact. 



b. Three qubits Take a special canonic wave func- 
tion with only four nonzero components: f/'oii; V'loii faw 
and ipooo — 1- The infinitesimal local transformations 
are 

5ip ijk = iPaicr^u'ipi'jk + iP^i^jj'^ij'k 

+ iP^((7l) kk 'faj k ' ■ (H7) 
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Let us show that Sipijk — implies P r K = for all 



K = ±,z, and r = 1,2,3. The 


corresponding 


system 


of equations is 






-i6i> 000 = -(Pg + P£ + P?) = 


o, 


(118) 


-iSipioa = P^ipiw + P^ipwi + 


p 3 - = , 


(119) 


—iSipoio — P^ipno + Pi~ipou + 




(120) 


— iSlbnm = Pnlbim -4- Pnthntt -+- 


P7 = 


(121) 


-i<fyon - (P 2 Z + Pi - Pg)ihu 


= o, 


(122) 


-^ioi = (Pf+P 3 *- p l)0ioi 


= o, 


(123) 


-z^ 110 = (P 3 2 +P 2 z -P 1 z )^iio 


-o, 


(124) 


-i^m = ^T^oii + P^i'wi + 




(125) 



Eqs. H119|l - I|121(l have the form of a homogeneous linear 
system of six equations for the real and imaginary parts of 
P, + . One can be convinced that the determinant of this 
system does not vanish for generic V'oii) ipwi, ipiw and 
this implies that the only solution is P+ = 0. Eqs. (|122|1 - 
(|124fl give in turn a homogeneous system for the three 
parameters P* = 0. Its determinant does not vanish for 
nonzero V'oiij "0101 > ipno and this implies P* = 0. Q.E.D. 

Note that this statement would not be correct for any 
pure state. The state having ^ooo as the only non- 
vanishing component is annihilated by two linearly inde- 
pendent generators (with Pf + Pf = and Pf + P 3 z = 0). 
For some other choice, there is only one trivially acting 
generator. But for a generic state there is none. 



c. n > 3 qubits The local transformations involve 
in parameters Pj\ Introduce, as before, complex P- = 
P* ± iPf ■ Take the wave function with O oo = 1 and 
n other nonzero components ^>oi...i> 0ioi...i, ■ ■ ■ , ipi...io- 
Consider 5ipox...i, etc. We obtain a system of n linear ho- 
mogeneous equations for n real variables Pf . The matrix 
of the system 




(126) 



(cf. Eqs. (|122ll - 1124() s ) has a non- vanishing determinant 
for n > 2. Indeed, M = A — 21, where A is the matrix 
with all components equal to 1 and / is the unit matrix. 
But A is a matrix of rank 1. It has n — 1 degenerate 
eigenvalues A = and one eigenvalue A = n. And A = 2 
is not an eigenvalue. It follows from this that the matrix 
M does not have zero eigenvalues and the only solution 
to the equation system is P z = 0. 

Consider now 5i/}if.\, where the set {k n } involves 2 ze- 
ros and n— 2 unities. We obtain a system of n(n— 1)/2 lin- 
ear homogeneous equations for n complex variables P^~ . 
Let us prove by induction that it has only zero solutions 
for generic 0oi...ii e ^c Consider first 5ip{k n }i with the 
unity at the leftmost position ({k n } — {1, {fc„_i}}). The 



equation system for P^ , . . . , P+ derived from this has 
the same form as the equation system <5^'{fc„_ 1 } in the 
case of n — 1 qubits. By the inductive assumption, it 
has only zero solutions. Knowing that P 2 ~ n = and 
considering, say, the equation 

<W>ooi...i = Piipio...i + -P 2 + 0oi...i = , 

we derive that in addition that P± = 0. 

This proves that no infinitesimal transformation leav- 
ing the generic n-qubit state exists for n > 2. 



d. Remark on the invariance of the canonic form 
Let us consider the canonic state suggested by Eq. @ and 
try to prove that this form is unique. Consider infinitesi- 
mal local transformations 1 + i J^i Pi a i + 1 Si Pi ~°i 01 
a canonic state vector and require that they do not bring 
about linear in cr^~ terms. The requirement means that 

#10.. .0 = ■ ■ ■ = 5ipo...oi = 

and implies 

P11P0...0 + ^ 2 + 0iio...o + • • • + P+1/V..01 = (127) 

P^ipia. . m + P 2 i'oio.. m + ■■■ + Pn^o...o = 0. 

This is a system of n complex equations for n complex 
parameters PJ = (Pj + ) , which we have encountered in 
Eq. I|65|) . It can be given in the form 

n 

Pr + Y J PtM id =Q, (128) 

i=l 

which with the allowance for Eq. I|61|l coincides with the 
right hand side of Eq. (|64|l . which we have encountered 
discussing the dynamic reduction of nilpotentials to the 
canonic form. There we have seen that all eigenvalues of 
Mij are never equal 1 for generic states in the canonic 
form, since it would contradict to the requirement of the 
maximum vacuum state population. Therefore the de- 
terminant of the system Eq. I|128|) does not vanish in the 
generic case, and hence the system Eq. <|127fl has only the 
trivial solution otj — 0. The canonic form of the generic 
state is thus unique and can only experience phase trans- 
formations, unless the phases are specified by additional 
requirements. 

B. Graph approach. 4-qubit si-classification 

In this Appendix we attribute a physical meaning to 
the coefficients of the tanglemeter by relating these to 
concurrence and 3-tangle. Based on these results, in the 
first section of this Appendix we construct graphs that 
illustrate the entanglement topology. In the second sec- 
tion, we summarize the s/-entanglement classes for four 
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qubits as they emerge from the dynamic equations for 
the tanglemeter. 



a. Graphical interpretation of the tanglemeter coef- 
ficients Starting with two qubits in a pure state, the 
tanglemeter involves only one coefficient, /?3, in terms of 
which we can express the von Neumann entropy S v n — 
— Tr[pA log pa] or the concurrence, by 

Cia = 2 v / det^ = 20 11 /{1 + /&) . (129) 

Since the tanglemeter in this case is nothing but the 
Schmidt decomposition, nothing new is introduced; a 
nonzero (3 coefficient implies the presence of bipartite en- 
tanglement among the two qubits. 

For three qubits in a pure state, both C and r can 
be expressed in terms of the amplitudes tp of the su- 
canonical state Eqs. Q291 I3l)fl . These expressions are 
shorter than for a non-canonical state vector, since they 
do not contain the amplitudes 0ioo, 0oiOi and 0ooi- 

For more than three qubits one can evaluate the con- 
currence. Since the expressions are rather complicated, 
we shall instead employ A. Peres' separability criterion of 
partial transposition |2q | : the qubits are not entangled iff 
the eigenvalues of the partially transposed reduced den- 
sity matrix of qubits 1 and 2 are non negative. The first 
case we considered is that in which only bilinear terms are 
present in the canonical state of four qubits. The eigen- 
values of the partial transposed reduced density matrix 
are found to satisfy the following relationships: 

K1K2 = |0oioi0iooi + 0ioio0ono| 2 - I ''/'ooiiV'oooo 1 2 , 
k x + k 2 =2Re (V>iooi0ioio + 0oioi0ono) , 

K3K4 = - |0oioi0iooi + 0ioio0ono| 2 

+ |0ooii0oooo 1 2 + I0OOH0HOO | 2 , 

«3 + k 4 = |0oooo| 2 + |0oon| 2 + |0noo| 2 (130) 

In the second case we consider, only trilinear terms are 
present, and we find 

kik 2 =, 100000 00m | 2 + |0oooo0ioii | 2 - I0111O0HO1 1 2 , 
ki + k 2 = |0oooo| 2 + |0ion| 2 + |0om| 2 , 

K3K4 = |0mo0iioi| 2 , 
k 3 + K4 = |0iiio | 2 + |0iioi | 2 • (131) 

The considerable simplification achieved for the for- 
mulas of familiar measures can help us to construct a 
topological picture of entanglement based on the invari- 
ant coefficients of the canonical form. This analysis also 
suggests an alternative derivation of the 3-tangle mea- 
sure. 

The rules that we prescribe for constructing the graphs 
to depict entanglement based on the coefficients of the 
tanglemeter are as follows: 



• Assign for each bilinear term a line connecting 
the qubits i and j. 

• For each trilinear term 0ijk, a surface on the plane 
confined by the k qubits. 

• For higher order terms, volumes, among the qubits 
involved. 

From Eqs. I|129[) . l|13U|l . one sees that the exis- 

tence of a bilinear term in the general case implies the 
presence of bipartite entanglement between the qubits i 
and j; we call this type of bipartite entanglement di- 
rect and represent it graphically as in Fig. El (a). Bi- 
partite entanglement is also present when the two qubits 
are indirectly connected by a line that passes through a 
third qubit. We can see from Eq. (|29|l that the concur- 
rence C12 is nonzero also when both 05,06 are present, 
Fig- (b). On the other hand, there is no indirect bipar- 
tite entanglement if the line connecting the two qubits 
involves more than two edges, as in Fig. El (c). In the 
general case, in which both direct and indirect contribu- 
tions are present, there are cancellation effects instead of 
addition, since the closed loops contribute to tripartite 
entanglement (see below). There are also cases in which 
bipartite entanglement is due to higher-order terms, as in 
Fig- El (b), where the terms /3mo, /?noi are present and 
the eigenvalues K\ and k 2 in Eq. I|131|l become negative. 
Another configuration of surfaces, like that in Fig. El (c), 
does not create bipartite entanglement. 
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(a) (b) (c) (d) 



FIG. 5: Bipartite entanglement among the qubits 1 and 2. 
(a) direct, (b) indirect, (c) configurations without entangle- 
ment, and (d) configurations that can result in cancellations 
of entanglement. 

Trilinear terms correspond to surfaces and by the 
Eq. (|29() we see that their presence in the canonical form 
of three qubits results in nonvanishing 3-tangle, as in 
Fig. El (a) . This is not the only configuration that per- 
mits tripartite entanglement; a loop consisting of three 



35 



lines is also a configuration with genuine tripartite en- 
tanglement, as in Eq. (|29|) . The inverse statement is not 
true for three qubits, since the term fa is not present in 
the expression of concurrence. 




<-i) (b> 



FIG. 6: (a) Tripartite entanglement in two different configura- 
tions, (b) The three classes of entanglement for three qubits. 
B stands for biseparable, W for the singular class containing 
the Werner state and G for the the general orbit. 

As an application, in Fig. [fJl (b) we represent the three 
classes of entanglement for three qubits diagrammatically 
way. It is important to note that if a diagrams is con- 
nected, i.e., there is a closed loop or surface, then this 
property cannot change under general local transforma- 
tions. 

We anticipate that the diagrammatic technique can be 
generalized to more qubits and improved with directed 
lines, surfaces and volumes that would aid in visualizing 
of cancellations or joint contributions. It appears that TV- 
partite entanglement can be attributed not only to the 
coefficients in front of the iV-th order monomials, but 
also to products of lower-order coefficients. 



b. Summarizing entanglement classes under SL{2, C) 
transformations In Sect. III.D.4, we have shown how do 
the s^-entanglement classes emerge from the considera- 
tion of the tanglemeter dynamics under controlled action 
of continuous local transformations when proper feed- 
back requirements are imposed. The general class has 
been identified along with four different singular classes, 
corresponding to one, two, three, and four vanishing 
eigenvectors of the determinant JHU- All these classes 
were characterized by three complex parameters. Here, 
we further refine this analysis and add some singular 
classes that depend on less parameters in order to com- 
pare this classification with that given in Ref. |24|. We 
illustrate the correspondence, when normalize the states 
of Ref. |24j to unit vacuum amplitude and compare these 
states with ones suggested by sZ-tanglemeters. 

Let us remind that the sZ-tanglemeter for four qubits 
we obtain in two steps. We first employ eight out of 



twelve complex parameters of local SL transformations 
and put the generic state depending initially on fifteen 
complex parameters to the form of Eq. (|27|) . which de- 
pends only on seven complex parameters. Then we chose 
the rest of the available transformations, four scaling op- 
erators e Bi(7i , in order to reduce the sZ-tanglemeter to the 
form of Eq. (|28fl. However, the second step not always 
possible to perform: when the tanglemeter is singular 
and one or more of the coefficients f3 equal zero, the sear- 
ings can simplify the nilpotential further, although not 
to the form of Eq. (J2SJ|. For example, when in Eq. lt77|) 
fa = 0, the tanglemeter can be set in the form 

fc = 4 CT 4 + + + CT 2 + 4) 

+ fa{cr+cr+ + CT+CT+) + o+o+a+(j+ , 

characterized by only two parameters. For fa = fad = 
it reads 

fc = °t a i + cr+ (7+ + fa(<jfai + <J^a+) + cr+a+(T+(T+ , 
and fa = fao = fa = results in 

f C = (T+a+ + (7+CT+ + (T+CT+ + (7+ (7+ (7+ (7+ . 

When the four-linear coefficient fas — and one or 
more of the quadratic coefficients are also zero, the sin- 
gular classes of states without fourpartite entanglement 
emerge: the sZ-tanglemeter of four-qubit W state 

fc = et a t + a i a 3 + a 2 a t^ 

belongs to one of these classes, and separable states with 
tanglemeters of the type 

f C = vjvj + + a 2 a X ! 

and similar, belong to other. 

Consider now several special cases of the singular class 
Eq. (|75|) , where all four of the eigenvalues of the determi- 
nant I|b9|) are zero. When, in addition, one or several of 
the trilinear coefficients equal zero, it becomes possible 
to rescale one or more of the bilinear coefficients to unity. 
Less general classes with sZ-tanglemeters like 

fc = vivi^i + oi<rt a i + ^t^i^i + a t°i 

fc = G i G 2 a i + G i G t a X + a l a 2 + <7 1~ £7 3 + P6 G 2 a 3 ) 

fc = <rt a £ a 't + G f G t + vi^t + a t a t 

emerge as a result. One reveals more singular classes, 
when besides of several zero cubic coefficients, two or 
more quadratic terms vanish. Two examples 

fc = G l G t a i + a i (T 2 rT 4 + °2 a £vi . 
fc = vtvivt + °x°i a X + ■. 

illustrate this case. 
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Different four-qubit classes emerging from this classifi- 
cation are presented in Table IWEl along with the results 
of Ref. [2J|. We include the general and main singular 
classes and omit a number of separable singular classes 
corresponding to product states. For the classes not sym- 
metric under cyclic permutation of qubit indexes, this 
permutations is implicit. We note that our classes result 
from consideration of dynamic evolution that implies a 
series of sequential infinitesimal local operations preserv- 
ing su-canonic form of the nilpotential. Therefore the 
situation, where some of the obtained classes turn out to 
be equivalent under a finite local sZ-transformation, yet 
cannot be excluded with certainty. Keeping this in mind, 



it is easy to see that G„. is identical to G a b c d class that is 
also suggested in Ref. £3 • The class L a b C2 corresponds to 
LG2 a , while L a2 b 2 coincides with LGl a - After applying 
SL transformations on qubits 3 and 4 of the L a b 3 class, 
the latter reduces to a form that is a singular case of the 
general class. The state L a4 can be set in the canonic 
form by flipping the second and third qubit, and then 
it is a special case of G c for f3§ = and /3j = i. Con- 
tinuing, L ao r>^„, coincides with St, Ln,»T with Sb, and 
Lo mI o 3fSl with S d - The singular class £o 5EB3 is of the S a 
form. Thus, we can conclude that the two classifications 
do overlap and complement each other. 
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C. Equation for local and two-body interaction 



where the coefficients 



We give here some more details on the derivation of 
the dynamic equation for the nilpotential /. Consider a 
single gate operation applied to qubits i and j. We cast 
the function / in the form foo+of fm+rf fw + ^t fn, 
and find 
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e~f = e-' 00 [1 - aff 01 - a+ f w - n+a+{hi - Joi/i 
c f = c /oo [j + a +f Q1 + a + fio + a + a +( fll + / 0l/l0 )] , 

(132) 



/oo = / - vJ - <?7 ~ i~ + (rja 



dot ~ 3 da+ 1 3 da+da+ ' 



, df + a 2 / 

JOl = ~ i~ - cr + 



3 da+da+ 



. df 



d 2 f 



dafdaj 



hi = 



d 2 f 



dajdaj 



(133) 



are independent of <jf , er^~ . Upon substituting these bearing in mind Eq. (|54|l , we obtain 
equations into Eq. Il- r >()[l for the Hamiltonian in Eq. Ij58(l . 



'M- = P-(T + + P-<J+ + P+-^L + P+^L - P a 

dt 1 1 + 1 dat 3 dat 1 1 



dat 3 3 dat 



Gi 




(134) 



Summing over z, j yields Eq. (|59|) . where the condition 
Gij = Gji is taken into account. 



D. Derivation of 3-tangle based on the tanglemeter 

Here we present a derivation of 3-tangle that is based 
directly on the transformations of amplitudes under lo- 



37 



General class 3 complex parameters 

Go / = f3 3 (a+a+ + a+a+) + f3 5 (a+a+ + a+a+) 

+P 6 (at at + at at) + a i at a z a t 
Singular 3D classes 3 complex parameters 

G b f = 3 (CTJ+CT+ + a+a+) + (3 B {afo+ + <T+(J+) + /3 6 (cr+cr+ + a+a+) 

+&1 at °t - a i °t °t + a i °i °t - a i a i °t 

+ 2 (#5/36 - PzPq + fefe) O+O+G+O+ 

+07 {at °i °i - °i °i at) + Pn {at at at - at at at) 

Gd f = at a 2 + a t a 3 + a 2 a t + °i a t + 0\<*t + a 2 a t 

+/?i4 {at at at - G t + a t at at - at) 

+/3i3 (crt<r£<rt + a t a t a t ~ a t a t a t - a t a t a t) 

{<yt a i ff t - a t a t a t - a t a t a i + at at at) 

G e f — at at at + at at at + a t a t a t + a 2 a f a t + 

03<rt °t + fla°2 a t + Psat &t 

Singular 2D classes 2 complex parameters 



LG2 a f = a+a+ + P 5 {at at + a+a+) + f3 6 (ata+ + a+a+) + ata+a+aj 
LG2 b f = atot + a+ + /3 5 (>+cr+ + a+cr+) + /3 6 (at°t + a£a£) 

LG2 C f = at at at + at at at + at at at + at at + feat^t + P^t^t 
Singular ID classes 1 complex parameters 



LGl a f = at at + at at + Pe{ata£ + a t a t) + at at at °"4 

LGlb f = at at at + at at at + at at + at at + Peat at 

Singular point classes no parameters 

S a f = at at at + at at at + at at 

Sb f = at at at + at at at + at at at 

S c f = a£ at at + at at at 

S d f = at at at 

S e f = at at + at at + at at + at at at at 

S f f = at at + at at + at at + at at at at 



TABLE I: Classification of four-qubit entanglement classes following from SL(2, C) transformation properties of the canonic 
form. 



cal operations. The explicit expression for 3-tangle has 
the rather simple form of Eq. (|30|l , when written for the 
su-canonical state. When considering SL transforma- 
tions, we look for a quantity that takes different values 
within the generic si orbit of three qubits, and vanishes 



outside this orbit. Performing SX-transformations, we 
can still set appropriate conditions so as to preserve the 
su-canonical form of states. 

We start with the s^-canonical state, i.e., is the maxi- 
mum entangled state 
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( 1 000) + I 111)) /\/2, written in terms of nilpotent vari- 
ables and normalizing to unit reference state amplitude: 

IV) = 4(1 + ^44)1000) . 

Then we apply to each qubit the general SL(2, C) trans- 
formation, which reads 



e A,«r eBi „* e C,„2 



A^B^Ci eC,i = 1,2,3. 



This transformation results in the most general admissi- 
ble su-canonical form of the su-canonical state, 



IV) = (Vooo+Voii44+V4 



10 u 2 



+ V4oi44 + Vin444)|ooo) 



in which the coefficients take the forms 



Vooo 
Von 

V4lO 



( 1 + Z ) C -Bi-B 2 -B 3 



V2 
1 + z 

Zy/2 
1 + z 

1 + z 

1 + 2z 



(135) 



z- 



•■V2 



CiC 3 e Sl+S3 - B2 
C 2 C3e S2+i33 - Bl 
C 1 C 2 C 3 e Sl + B2 + B3 



Here we have employed a complex variable z that nat- 
urally emerges from the requirements of su-canonical 
form, having the relations z = Aie 2Bl C\ = A2e 1B ' 2 Ci = 
A 3 e 2B3 C 3 and (1 + z) 2 C 2 C 3 = -z 2 . One sees that the 
amplitudes of the sw-canonical state depend on seven 
real parameters, since in addition to five parameters of 
the tanglemeter, we now allow for two more parameters 
characterizing the vacuum state amplitude. 

Being in the general orbit, states of less general orbits 
can be reached if irreversible transformations are per- 
formed [35| corresponding to the limit |z| — > oo. At this 
point we want to construct a polynomial measure on the 
amplitudes tp', such a way that it is a function of z and 
it vanishes outside the general sZ-orbit. Direct inspection 
shows that we can construct a function that depends only 
on the parameter z 



C(*) 



Vll 1^000 



= (l + 2z) 2 
Vooo^iio^oiiVioi z(l + z) 



(136) 



and is independent of the normalization,. The function 
£(z) shown in Fig. [7] can take any value, but —4 that 



corresponds to the limit |z| — > oo. Therefore, the poly- 
nomial 



V'mV'ooo + ^oooViio^oiiV'ioi 



will be identically zero outside the general orbit, and is 
otherwise nothing but the 3-tangle expressed in terms 




FIG. 7: Function \((z)\ as in Eq. JT^. In the limit \z\ -» oo, 

|C(*)|-4 



of the coefficients of the canonical nilpotent polynomial, 
Eq. CP. 

This procedure for constructing a measure can be in 
principle extended to more qubits. However, the exten- 
sion will require some care. The general sZ-orbit of four 
qubits is characterized by six parameters parameters, in 
contrast to the three-qubits case. Consequently, irre- 
versible transformations may connect different sZ-orbits 
that both contain iV-partite entanglement. Therefore the 
desired limits need to be clearly specified. 

Finally, we would like to mention that Theorem 3 in 
Ref. [3^ suggests another useful application of the £(z) 
function. For an arbitrary 3-qubit state expressed in 
canonical form, one can calculate the value of £(z) by di- 
rect substitution of the numerical values of amplitudes to 
Eq. 113bj) , and then solve a binomial equation to find the 
root z — z r . One can choose the coefficients Ai, Bi, and 
Ci such that they satisfy Aie 2Bi Ci = z r . In this way a 
determinant- 1 transformation that brings the maximum 
entangled state to the chosen state can be identified ex- 
plicitly. The inverse transformation can be used for an 
optimal filtering procedure called purification, a proba- 
bilistic procedure that transforms a state of the general 
orbit to the one with the maximal entanglement, i.e. to 
the GHZ state. 
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